A longa distância (r > 29) predominam as forças de 
dispersão de London, polarização ou electrostáticas — 
atractivas. 


Não existe nenhuma expressão analitica que traduza, 
duma forma univoca, o potencial intermolecular, pelo que 
têm sido desenvolvidos vários modelos matemáticos que 
pretendem ser uma aproximação fisica ao potencial inter- 
molecular. 


Dos modelos desenvolvidos foram aplicados, no pre- 
sente trabalho, os modelos de Lennard-Jones e Kihara [9] 
dados respectivamente pelas equações 5 e 6, e o modelo 
de Dymond-Alder [10] que sendo um potencial numérico 
não tem expressão analítica. 


&(r) ( ai ME 

E od , si 
ME (ey Tom 

E À r r 


r=(0—d) V2. (6b) 


No potencial de Kihara d é o diâmetro do cerne 
rigido e a sua importância costuma ser expressa pelo 
parâmetro Y = d/(0—.d). 


3.2. Escolha dos parâmetros 
3.2.1. Método das isotérmicas 


O cálculo dos parâmetros E e q pode ser feito se- 
guindo o método das isotérmicas [9]. Este método ba- 
seia-se na utilização da equação (1) na forma 


5 (x mkT)!? 
= — — = (7) 
16 x 02 m(T) 


o 3» (1º) 
E (7º) 


em que Tº =kT/£ é a temperatura reduzida e no lado 
direito da equação (7) utilizam-se valores experimentais 
de m(T) para cada gás. Para se iniciar o cálculo torna-se 
necessário usar uma primeira estimativa de q para o gás 
em causa que foi escolhido comparando os valores de q 
obtidos por outras propriedades e também pela viscosi- 
dade, para cada tipo de potencial [9]. 

Admite-se que a (T*) = 1, e para cada par de pon- 


tos n(T) obtém-se Q (2.2)* (T*). Pela consulta dos inte- 
grais de colisão para cada potencial, obtém-se o valor 
de T* e consequentemente o valor de £. Ficamos com 
um par de valores (£, 0) a uma dada temperatura. Uma 
vez conhecido T* podemos reciclar, obter mn (T*), a razão 


(7) e obter um novo valor de E correspondente a q. 
Escolhem-se seguidamente outros valores de o e obtém- 
-se uma recta para cada temperatura, quando se repre- 
senta E em função de q. O processo é repetido para 
outras temperaturas obtendo-se um feixe de rectas, que 
se cruzariam num ponto único se o potencial utilizado 
fosse o potencial real entre as moléculas. Como na 
realidade se trata de um modelo diferente do potencial 
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real, as rectas cruzam-se numa zona do diagrama (E, 0) 
esboçado na Fig. 2. Dum modo geral, os dados de baixa 
temperatura conduzem a valores baixos de E e altos de q 
enquanto que os dados de alta temperatura conduzem 
a valores altos de € e baixos de 0. 

A escolha dos valores de (£, 9) recai sobre um par 
de valores médios, os mais prováveis da zona referida. 


Elk 


Emédio 


Omédio 


FIG. 2 


Cálculo dos parâmetros dos potenciais a partir dos dados 
de viscosidade T, <T; <T; <T, 


Este processo foi utilizado para os potenciais de 
Lennard-lones 6-12 e de Dymond-Alder, no primeiro caso 
para o crípton e azoto e no segundo para o azoto e 
metano, encontrando-se os valores obtidos na Tabela |. 
O valor dos integrais de colisão foram obtidos para os 
potenciais de Lennard-lones e Dymond-Alder da referên- 
cia [9], e para o potencial de Kihara da referência [14]. 


3.2.2. Método dos estados correspondentes 
Outra forma de obter os parâmetros é a baseada no 
princípio dos estados correspondentes, em que a razão 


entre os parâmetros é função das razões entre as pro- 
priedades críticas, isto é, Te, Vc e Pc [11]. 


(: (: ) (To); 
=) “Era (8) 
k : k Ar (Todar 


(Vol; 


(bo); = (Doar - (9) 


Char 


em que 


bb=— 7 No (10) 


vim 


Os parâmetros de potencial do árgon foram tirados 
de [9] para os diversos potenciais e os valores utilizados 
estio presentes na Tabela |. As constantes críticas dos 
diversos gases foram obtidas em (12). Para o potencial 
de Kihara foi admitido Yk, = 0,1 e Ycm a a 0,3, valores 


das entradas das tabelas dos integrais de colisão do 
potencial de Kihara [13], que mais se aproximavam dos 
parâmetros obtidos por TEE et al [14]. 
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4. RESULTADOS OBTIDOS 
4.1. Gases puros 


Após terem sido obtidos os parâmetros moleculares 
de cada gás para os diferentes gases, calcularam-se os 
valores da viscosidade teórica para o Kr, N; e CH,, bem 
como da condutibilidade térmica e do coeficiente de 
autodifusão do Kr e da condutibilidade térmica do me- 
tano. 

Os valores experimentais da viscosidade dos gases 
foram tirados da análise de MAITLAND e SMITH [15] 
não sendo valores experimentais «stricto sensu- mas 
valores obtidos dum ajuste de minimos quadráticos dos 
resultados experimentais existentes até à data e cobrem 
os seguintes intervalos de temperaturas: 


máximo erro em relação aos 
valores experimentais: 


Kr 140 K <T< 2000 K + 1,0% 
N; 100 K < T< 2000 K + 1,5% 
CH, 100K<T< 1000 K + 13% 


A razão t, varia pouco com a temperatura e com o 
tipo de potencial e segundo KESTIN [16] pode ser dada 
por: 
+ 
 (T*) = a = 1 + 0,0042 [1 — exp(0,33 (1 =: "98 (12) 
n 


Para a gama de temperatura utilizada a função & 
(T*) não se afasta da unidade mais de 0,3%, pelo que 
não foi utilizada nos cálculos, pois os resultados expe- 
rimentais utilizados para comparação têm um erro esti- 
mado em 1%. 

No caso do metano, dado tratar-se de uma molécula 
poliatômica, torna-se necessário aplicar uma equação 
para a condutibilidade térmica que entre em linha de 
conta com a possibilidade de existência de colisões ine- 
lásticas. A primeira tentativa de cálculo da condutibili- 
dade térmica dum gás poliatômico deve-se a EUCKEN 
[17], tendo a expressão que este autor deduziu sido 
modificada por UBBELOHDE [18] admitindo que a energia 
dos graus de liberdade rotacionais e vibracionais é per- 
mutada segundo um processo difusional. A expressão 


TABELA | 


Parâmetros moleculares do Ar, Kr, N: e CH, determinados pelo método das isotérmicas a partir de dados de 
viscosidade e a partir do principio de estados correspondentes, para os potenciais de Lennard-lones 6-12, Kihara 
e Dymond-Alder 


LENNARD-JONES 


Dados /T) P. E. C | 


ElkiK cinm ElkIK ginm elkiK 


o —T— 


Ár 120,72 | 0,3434 ao = 143,26010) 
Kr | 205,0 0,350 — ue 199,0 

N: 108,4 0,360 E | —— = 
CH. — — | 152,5 | 0,377 | 232,2014) 


819 | KIHARA 


0,351 414) 0,2848) [1740 


DYMOND-ALDER 


; Es. Dados m/T) P.E. C 

ginm | Y EkiK | sinm EkiK cinm 
0,3344009 | g11109 | 138,20 | c 328010) — — 
0,357 0,1 — — 1920 | 0,350 


“| > [mo [os ER 


a — | ST —— e À O SD 


0,360 174,7 0,360 


NOTAS: 1. Os números entre parêntesis referem-se à bibliografia donde foram tirados os valores para o caso do árgon, tomado como 


referência no método do P. E. EC. 


2. Para o metano os parâmetros obtidos pelo P. E. C. para o potencial de Kihara eram bastante diferentes dos tabelados('9) 


pelo que utilizámos os valores de TEE et al (1), 


Os valores experimentais da condutibilidade térmica 
do cripton foram obtidos das tabelas de VARGAFTIK [20]. 

Os valores teóricos para a condutibilidade térmica 
do cripton foram obtidos a partir das equações 1 e 2 
que podem ser combinadas na equação (11): 


15 R f 


RE o Dr 1 
4 M f 82 


3 
na medida em que para um gás monoatómico Cy = — RR. 
2 


162 


proposta por estes autores, válida normalmente para altas 
temperaturas é: 


( AM * R 
=1,77— + 1,32 13 
E U V Pa 


L 


MASON e MONCHICK [8] propuseram uma teoria 
mais correcta, mas a dificuldade em obter valores do 
parâmetro Z,y levou-nos a aplicar numa primeira fase, 
a equação de EUCKEN modificada [18], equação (13), e 
posteriormente a teoria de MASON e MONCHICK. 
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A equação de MASON e MONCHICK [8] introduz o 
parâmetro Z,y que é o número de colisões necessárias 
para haver permuta da energia rotacional e translacional, 
e traduz-se por: 


“AM E dO Sr 2,184 Cro 
de SEE E - (14) 
n Cy /MM NCy /uy Tt Cy Erik 


em que Cy é a contribuição rotacional para o calor 


específico do gás Tú para o metano). 


Na ausência de valores experimentais de Z,, ele 
funciona como um parâmetro ajustável. 


Os valores experimentais da condutibilidade térmica 
do metano foram obtidos de (19) e os do coeficiente de 
autodifusão de (20). 


Os resultados obtidos para os- diferentes modelos 
de potencial e para as diversas propriedades de trans- 
porte foram representados na forma de percentagem de 
afastamento do valor teórico do valor experimental, 
isto é: 


X teórico — X experimental 
a erro = — e 


? — x 100% (15) 
X experimental 


e apresentam-se nas figuras 3, 4 e 5 e nas tabelas !l e Ill. 


CRIPTON 
VISCOSIDADE 


CONDUTIBILIDADE TÉRMICA 


FIG. 3 


Viscosidade e condutibilidade térmica do cripton segundo o poter- 
cial de Lennard-lsnes — + . de Dymond-Álder — 6 . e de 
ronera — & . comparação com os resultados experimentais 
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AZOTO 
VISCOSIDADE 


500 1000 1500 T/K 2000 


FIG. 4 
Viscosidade do azoto segundo o potencial de Lennard-lones — + 
e de Dymond-Alder — O . Comparação com os resultados expe- 
rimentais 
METANO 


VISCOSIDADE 


10 õ 
E o 
o 
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FIG. 5 


Viscosidade do metano segundo os seguintes potenciais 

O — Lennard-lones, parâmetros citados na Tabela | (P. E CJ) 
6 — Lennard-Jones, outros parâmetros (P. E C) 

“e — Dymond-Alder, método das isolérmicas 

E — Dymond-Alder, (P. E. C.) 

e — Kihara, parâmetros da referência (14) 

O —Kihara, (PE. C) 
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TABELA Il 


Condutibilidade térmica do CH, gasoso segundo o modelo de Lennard-lones (6-12) e comparação com os resultados 
experimentais [20], para p=1 atm. Os cálculos foram efectuados pela equação de UBBELHODE e de MASON e 
MONCHICK, com Za = 1,06 


O OD) 


NmwWm-1K-] (20) o erro 
TK Te - À exp imWme IK-1 T 
UBBELHODE MASON/MONCHICK eq(13) eq(14) 
335,6 2,20 46,18 38,31 38,31 20,5 0,0 
371,3 2,43 52,59 44,05 44,00 19,5 0,1 
403,8 2,65 58,66 49,53 51,12 14,7 = 3,1 
435,6 2,86 64,64 54,97 53,59 20,6 + 2,6 


TABELA II 


Coeficiente de autodifusão do CH, gasoso segundo o modelo de Lennard-lones (56-12) 
e comparação com os resultados experimentais [21], para p=1 atm 


TIK 


273.16 


298,16 


105D xp/mês-] Sa erro 
2,06 | —9,7 


2,40 —8,8 


4.2. Mistura Argon-Cripton 


O coeficiente de difusão mútua do argon em cripton 
foi calculado através da expressão de Chapman-Ens- 
kog [3]: 


Dya/em? "= 
VIM, + MM, M; 
= 26280104 — em PTOS  ftB) 
2 OU,DIT* Dt 12 
(p/atm) o pc do) 


em que M, é a massa molecular do constituinte |, e a 
expressão já está preparada para utilizar o valor da 
pressão p em atmosferas. Para podermos aplicar a equa- 
ção (15) necessitamos dos valores das constantes de 
força Cj e €2/k para a interacção árgon-cripton. Num 
primeiro passo testaram-se as regras de combinação de 
Lorentz (eq. (16)) e Berthelot (eq. (17)), que vêm dadas 
por: 


1 
Cj = E (0, + 03) (16) 


Elk = Vey)k . Es/k (17) 


utilizando os valores dos parâmetros do árgon e do 
crípton constantes da Tabela |, para o caso do modelo 
de Lennard-Jones, tendo-se obtido ems/)k = 157,3K e 
4 = 0,347 nm. Os valores obtidos podem comparar-se 
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com os parâmetros deduzidos por KESTIN et al [21] para 
o potencial de Lennard-Jones a partir do ajuste aos valo- 
res experimentais dos coeficientes de difusão, respecti- 
vamente 171,67 K e 0,3408 nm, o que será feito na secção 
seguinte. Para ambos os pares de parâmetros foi cal- 
culado o valor do coeficiente de autodifusão no intervalo 
de temperaturas 200-1000 K que se comparou com os 
valores experimentais acessíveis [20]. Os resultados obti- 
dos estão representados na figura 6. 


5. DISCUSSÃO DOS RESULTADOS OBTIDOS 


A teoria de CHAPMAN-ENSKOG é, como se refere 
na introdução, válida para gases monoatómicos, sendo 
uma teoria exacta caso seja conhecido o verdadeiro po- 
tencial intermolecular, o que não acontece na realidade. 

Para a cripton, uma vez determinados um conjunto de 
parâmetros E/k e q, consistente com um dado modelo de 
potencial, a partir de uma propriedade, as outras proprie- 
dades de transporte podem ser calculadas com uma 
exactidão razoável. O afastamento dos valores calculados 
em relação aos valores experimentais traduzem a imper- 
feição do modelo de potencial e as suas limitações em 
traduzir o verdadeiro potencial intermolecular. 

Assim, dos três modelos de potencial utilizados para 
o cripton, apenas para o modelo de Lennard-lones se 
determinaram os parâmetros duma forma absoluta (mé- 
todo das rectas isotérmicas), verificando-se que a visco- 
sidade é sobestimada em 2% em todo o intervalo de 
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200 300 400 500 500 700 800 


900 1000 T/K 


FIG. 6 


Coeficiente de difusão mútua da mistura argon-cripton. 
o — Lennard-jones com regras de combinação 
* — Lennard-lones, parâmetros de Kestin et al (21) 


a — Valores experimentais (22) 
e — Valores experimentais (23) 


temperaturas utilizado, enquanto que a exactidão dos 
valores experimentais utilizados é de + 1%. À conduti- 
bilidade térmica é reproduzida pelo modelo de Lennard- 
-Jones a + 1% para 273< T< 800K divergindo a bai- 
xas e altas temperaturas. Os potenciais de Kihara e 
Dymond-Alder subestimam os valores da codutibili- 
dade térmica, aumentando o afastamento com a tempe- 
ratura. Para os potenciais de Dymond-Alder e Kihara 
foram utilizados parâmetros relativos calculados pelo 
princípio dos estados correspondentes, sendo a concor- 
dância razoável (—1,5 %) para o potencial Dymond-Alder 
no caso da viscosidade (só determinado com valores de 
propriedades do árgon), embora para o potencial de 
núcleo rígido de Kihara o desvio aumente com a tempe- 
ratura. 

O azoto é uma molécula diatômica com uma energia 
rotacional elevada pelo que se comporta praticamente 
como esférica (valores absolutos pequenos do 2.º coefi- 
ciente de virial, por exemplo). Sendo assim, não é de 
espantar o excelente acordo obtido para o potencial de 
Lennard-lones, com os parâmetros consistentes com a 
propriedade para a qual foram ajustados. Salienta-se, 
contudo, a sua capacidade em reproduzir os dados expe- 
rimentais da viscosidade a menos de 0,3 % num intervalo 
da temperaturas enorme (300 a 2000K). No caso do 
potencial de Dymond-Alder este acordo já não se verifica 
para o conjunto de parâmetros utilizados, calculados pelo 
princípio de estados correspondentes. 
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Para o metano, a molécula globular, de cume simé- 
trico, mais simples, qualquer dos potenciais moleculares 
utilizados reproduz os dados da viscosidade a + 2% 
no intervalo de temperaturas estudado (300 a 1000K), 
divergindo grandemente para baixas temperaturas onde 
a importância das forças repulsivas se torna preponde- 
rante. À única excepção é o potencial de Kihara com 
parâmetros moleculares determinados pelo princípio de 
estados correspondentes em que o afastamento entre 
os valores calculados e os experimentais é grande (6 %). 

No que se refere à condutibilidade térmica e ao 
coeficiente de autodifusão (Tabelas Il e III) apenas foi 
testado o modelo de Lennard-Jones; a equação de UBBE- 
LOHDE [19] sobrestima em cerca de 20% o valor da 
condutibilidade térmica do metano e subestima em 10% 
o coeficiente de autodifusão o que pode ser devido ao 
facto de os parâmetros do metano terem sido calculados 
pelo princípio dos estados correspondentes, facto aliás 
suportado pelo cálculo do coeficiente de autodifusão 
(10% de desvio). 

No entanto, a equação de EUCKEN modificada 
sobrestima normalmente os valores da condutibilidade 
térmica de gases poliatómicos simples, como por exem- 
plo para o CF, [9]. A utilização da equação de MASON 
e MONCHICK, com a aproximação p D/ n— 1,32 [9], 
com o valor de Z,w = 1,06, reproduz bastante bem a 
condutibildade térmica do metano no ntervalo de tempe- 
raturas estudado. 
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No caso da mistura árgon-cripton, como se pode 
verifcar na figura 6, o acordo entre os valores experi- 
mentais e os calculados pelo modelo de Lennard-Jones 
com os parâmetros da mistura calculados segundo as 
regras de Berthelot e Lorentz, na zona em que existem 
dados experimentais, é excelente. Para temperaturas mais 
elevadas não existem dados experimentais acessíveis e 
só se pode comparar os resultados obtidos com os 
resultados teóricos para o modelo de Lennard-Jones obti- 
dos por KESTIN et al [21], com cálculo dos parâmetros 
E ec a partir da viscosidade da mistura. O afastamento 
entre os dois conjuntos de resultados teóricos aumenta 
com a temperatura, sendo de cerca de 5% à tempera- 
tura de 900 K. 

Podemos assim concluir que a utilização de parâme- 
tros moleculares determinados a partir de dados de vis- 
cosidade, independentemente do modelo de potencial 
intermolecular utilizado, é um processo rápido e expe- 
dito de previsão de propriedades de transporte dos gases 
à pressão normal, divergindo no entanto os cálculos para 
as temperaturas baixas (T< 250K). A utilização dos 
parâmetros determinados pelo método das isotérmicas 
conduz a resultados bastante razoáveis para todos os 
potenciais o que atesta a consistência interna do método. 


Os parâmetros determinados a partir das relações do 
princípio de estados correspondentes dão resultados pio- 
res, devido à simplicidade das relações utilizadas para 
a sua estimativa. 


Posteriormente à realização prática deste trabalho, 
chegaram ao conhecimento dos autores valores dos pa- 
râmetros moleculares determinados a partir de técnicas 
de inversão [25], [26] e para o caso dos gases raros, 
a partir dos trabalhos de BARKER [27], calculados a 
partir de dados experimentais de várias propriedades, 
nos estados sólidos, líquido e gasoso. Os valores obti- 
dos duma forma absoluta por estes autores estão repre- 
sentados na Tabela IV e pode-se verificar o acordo 
bastante bom entre os valores de o para as diferentes 
moléculas correspondentes ao potencial de Lennard-Jjo- 
nes, incluindo a mistura árgon-cripton, embora para a 
dimensão do fosso de potencial o afastamento seja ligei- 
ramente maior. 

Dada a facilidade de cálculo de parâmetros inerentes 
ao método das rectas isotérmicas, recomendamos a sua 
utilização para estimativa rápida de parâmetros, para 
todos os problemas em que se requer uma precisão ele- 
vada. 


TABELA IV 


Parâmetros moleculares obtidos por diferentes métodos e sua comparação com os calculados no presente trabalho 


LENNARD-JONES 


DYMOND-ALDER 


MÉTODO DE BARKER (2) TÉCNICAS DE INVERSÃO 


6-12 (24-26) 
ElkiK g'nm 
0,357 199,2 (b) 0,359 (b) 
nas 217 (b) 0,356 (b) 
157,3 (a) 0,347 (a) 
cm 169,4 34 1 
171,7 (d) | 0,348 (d) 69,4 (a) 0,346 (a) 65 (c) 0,349 (c) 


(a) — regras de Berthelot e Lorentz 

(b) — SMITH et al (24-25) 

(c) — G. C. MAITLAND, W. A. WAKEHAM (26) 
(d) — KESTIN et al (21) 
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NOTICIÁRIO 


(Continuação da pág. 144) 


FEIRA DE LEIPZIG — REPÚBLICA DEMOCRÁTICA ALEMÃ — 9/3-15/3 — 1980 


Contactos a nivel mundial para um comércio pleno de êxito... 


...São-vos oferecidos pela Feira de Leipzig. 


Leipzig cria para si as melhores possibilidades para uma troca de experiências e comparação de resultados, 
especialmente nos campos da construção de máquinas e instalações, da metalurgia, da electrotécnica/electrónica, da 


construção civil e, ainda, no campo dos bens de consumo. 


Uma informação actual a nível internacional através de simpósios, colóquios e conferências especializadas. 


Esperamos por si em Leipzig! 


Informações adicionais e Cartões de Visitante poderão ser obtidos na 
Embaixada da República Democrática Alemã — Departamento Comercial 


Rua Barão de Sabrosa, 213 — Lisboa 


bem como em qualquer posto de fronteira da Repúblca Democrática Alemã. 


(Informação do «Expressos, de 9-2-1980) 


ACÇÕES DE FORMAÇÃO/1980 — CALENDÁRIO DAS ACÇÕES INTEREMPRESAS 


Cursos em Lisboa: 


Aperfeiçoamento de Preparadores e Encarregados — Inicio em 10 de Março. Duração: 10 dias. 
Gestão da Produção — Início em 7 de Abril. Duração: 5 dias. 

Téc. Lubrifcação Orgãos de Máquinas — Início em 14 de Abril. Duração: 3 dias. 
Organização de Armazéns e Gestão de Stocks — Início em 5 de Maio. Duração: 5 dias. 


Cursos no Porto: 


Organização da Conservação — Início em 25 de Fevereiro. Duração: 5 dias. 
Organização da Lubrificação — Início em 3 de Março — Duração: 3 dias. 
Informática Operacional na Área Técnica — Início em 10 de Março. Duração: 5 dias. 


DEPARTAMENTO DE FORMAÇÃO — Rua Bernardo Lima, 48-1.º — LISBOA 


Telefones: 55 46 95/57 93 77. 


TÉCNICA 455 


(Informação do «Expresso», de 9-2-1980) 
(Continua na pág. 188) 


167 


OLHOd +9+089 Jelel Obg 'elbaly olesseg op 'H | [| > I X HF | N a 
| As 
os sumo de MN USANINIAOA 9NOLA É - 


Do UUI/[2D% GS LO = dy - ODINHIL OLNIWNVTOS! S9NOLA Oy HOlHIdNS < (0V5d9VHL) VIONILSISIH 
LD S'7 - VHNSSIdASI SOrINZY 3 DNOLA VIHYNIATV - OVÔVZMILN 
SO 3 +'0 - aQVaISNIG SV100 € 


SVOINHIL SvHHOS é 
Do ULI/[29Y GSL'O = dy - ODIINHIL OLNINVTOSI 


DNOLA SVIHVYNIATV Sy TYNDI - OLNINVAVOV cUO/0Y 000'2L - IJAVAIDILSVTII OINNAOW 
(OZL Y) (otz V) 
SvHOH z Y HOIHIANS - 0904 OV VIONILSISIH SvHOH + V HOlHIdNS - 0904 OV VIODNILSISIH 
Do UUI//29% O2'0 = dy - ODINHIL OLNINVTOSI UO/DY 0€ - (OVSSIHANOD) VIONILSISIH 
TVNOIDIaVHL Ov HOIHIANI %O0E V OL - OSId Wu | + - TVNOISNINIA VIONVHITOL 
2LU/BY 00S - VINIXVIA VOHVO S 0 - JAVaISNIA 

SOHLIIN 00 6 - OWIXVIN OVA S09018 € 


OdvdlsoO!Iv OLNINIAVd é 


TECNICA IV 


NUMERO 455 DEZEMBRO 1979 


ANO LIV 


VOLUME XLI 


Relatório sobre a disciplina 


SISTEMAS DE CONTROLO E DE COMUNICAÇÃO 


do 5.º ano do Curso de Engenharia Electrotécnica 


do Instituto Superior Técnico !” 
Janeiro, 1978 


SUMÁRIO 


Estuda-se neste relatório a disciplina Sistemas de 
Controlo e de Comunicação do grupo «Análise e Teoria 
dos Sistemas» do curso de Engenharia Electrotécnica. 
Começa-se por reflectir sobre o âmbito da Teoria dos Sis- 
temas, especificando para o grupo de disciplinas as maté- 
rias afins e o espirito prevalecente Percorre-se depois 
sumariamente aspectos do ensino da Teoria dos Siste- 
mas no curso de Engenharia Electrotécnica no IST, 
nomeadamente no seu Departamento de Teoria dos Sis- 
temas, Finalmente discute-se a natureza da disciplina de 
Sistemas de Controlo e de Comunicação concretizando-se 
o respectivo programa, o conteúdo, métodos de ensino e 
aspectos bibliográficos. 


1. Teoria dos Sistemas: Uma Classificação 


Nos anos 40 Norbert Wiener propôs [1] que se 
baptisasse de Cibernética uma nova ciência, síntese de 
diferentes disciplinas como o Controlo e a Comunicação 
em diversos campos de aplicação, A figura 1 representa 
em diagrama um possivel retrato na época da génese da 
Cibernética. A Cibernética é a tomada de consciência 
da unidade prevalecente nos Computadores, no Controlo, 
na Comunicação e vai construir-se basicamente sobre a 
Electrónica que, constituindo o suporte material tecnoló- 
gico dos sistemas noutros campos, lhes impõe uma diná- 
mica acelerada de transformação que, rapidamente, de- 
sactualiza a fotografia anterior, Sobre uma perspectiva 
desta evolução ver a Parte | de [2]. 


Com toda a subjectividade inerente a qualquer sis. 
tema de classificação, propõe-se, como actualização nos 
anos 70 da mesma preocupação de Norbert Wiener, o 


JOSE MANUEL FONSECA DE MOURA 


Universidade Técnica de Lisboa (1.5.T.) 
Centro de Análise e Processamento de Sinais (I.N.1.C€.) 


ABSTRACT 


The paper presents Control and Communication Sys- 
tems, a subject of the group «Analysis and Theory of 
Systems» of the Electrical Engineering program. If starts 
out by reflecting on the nature of Systems Theory itself, 
specifying for the above group of subjects its contents 
and prevailing attitude. The report then summarizes as- 
pects of the teaching of Systems Theory in the Electrical 
Engineering Course at IST, namely at its Systems Theory 
Department. Finally, the paper discusses Control and 
Communication Systems presenting its programm, its 
contents, the pedagogical methods used and the required 
bibliography. 


esquema que se desdobra pelas figuras 2 a 8 O con- 
traste entre os dois diagramas reflecte a interpenetração, 
a expansão e o ressurgimento decisivo de novas discipli- 
nas e ciências. 


| Cibernética | 


4 


Controlo Comunicação 
Computadores 
e ae A 
Electrônica 
FIG. 1 


Cibernética de Norbert Wiener 


(!) De acordo com o artigo 9.º, parágrafo 1, alinea a) do decreto n.º 301/72 de 14 de Agosto. 


* Manuscrito recebido para publicação em 7-6-78. 
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Como nó central na figura 2 aparece a Teoria dos 
Sistemas. Deste emergem três ramificações que isolamos 
em qualquer aplicação: 


|) Compreensão dos fenómenos físicos subjacentes 
à experiência e sua teorização. 


ii) Implementação e seus aspectos tecnológicos, 
onde se singulariza a Electrónica como principal 
e mais fecundo meio de implementação. 


iii) Processamento, entendido como conjunto de clas- 
ses formais de operações comuns a diferentes 
contextos físicos e implementações. Inclui o esta- 
belecimento e análise do modelo formal e a sin- 
tese das acções de intervenção, Ou, dito de outra 
maneira, corresponde ao estudo e ao projecto 
formal dos sistemas e sinais que interactuam na 
experiência. 


Teoria dos Fenômenos 


Teoria 


dos 


Sistemas 


Processamento Tecnologia 


FIG. 2 


Teoria dos Sistemas: Suas Ramificações 


Às duas primeiras categorias são inerentes os con- 
ceitos de energia e de matéria. Na terceira é a informa- 
ção, e o seu fluxo, formalizados por Shannon em 1948 [3]. 
Às duas primeiras preocupações são por isso o objecto 
privilegiado de disciplinas decorrentes da Fisica, da Qui- 
mica, enquanto o processamento tem paladar de natureza 
formal, ou seja, matemático, 

Um Exemplo: Antes de prosseguir, esclarece-se o 
diagrama proposto com uma breve discussão centrada no 
contexto da Comunicação. Este vocábulo aparece tradi- 
cionalmente associado à transmissão a distância (teleco- 
municação) de mensagens escritas, orais ou visuais, em 
que o homem é num extremo a fonte e no outro o des- 
tinatário. Nesta acepção a expressão é associada aos 
grandes sistemas de comunicação historicamente desen- 
volvidos — o telégrafo, o telefone e a telefonia (rádio), a 
televisão. 

A figura 3 ilustra em esquema de blocos um sistema 
de Comunicação, O mesmo meio fisico, p.e. na rádio e 
na televisão, transmite mensagens de diversa natureza. 
O mesmo tipo de mensagens, p.e, no radar e no sonar 
pode ser propagado por suportes fisicos diferentes. 
Nestes dois últimos exemplos o canal de transmissão, 
bloco central da figura 3 a), é o espaço livre ou o canal 
sonoro do oceano, respectivamente Em qualquer deles 
se observa o fenómeno da propagação, Mas, ao radar, 
interessam os fenómenos electromagnéticos, objecto de 
estudo do Electromagnetismo, enquanto ao sonar as ondas 
e a sua propagação são de outra natureza, do dominio 


170 


da Acústica. A especificidade dos fenómenos fisicos 
determina o tipo de interface (ver a figura 3-b) entre os 
blocos iniciais do emissor e o canal, pelo que os hidro- 
fones dum sonar se distinguem inequivocamente dos di- 
polos duma parabólica, 

Nos exemplos anteriores é possivel isolar operações 
que abstraiem da natureza da mensagem, da sua fonte, 


| Processador [ | InterFace 


Emissor 


tb) 
FIG. 3 


a) Sistema de Comunicação; b) Emissor 


do suporte físico do meio de transmissão, do destinatário. 
É o processamento que relaciona o conteúdo informativo 
com a dinâmica dos fenómenos, i.e. com as variações 
no tempo dos fenómenos observáveis. À natureza (fisica) 
do problema implica com os factores de ordem técnica 
de implementação. Mas, a determinado nível de abstrac- 
ção, há uma sucessão de operações conceptualmente 
equivalente, independentes da natureza do canal, O pro- 
cessamento surge, a este nivel, desvinculado do meio 
fisico subjacente. 

à universalidade no exemplo anterior é bem mais 
lata e ultrapassa o âmbito de diferentes aplicações no 
mesmo campo, De facto, o processamento, tal como 
entendido acima, pode ser isolado como fase intermédia 
em aplicações em campos totalmente diversos, Passa-se 
da Comunicação (na acepção tradicional) para a Econo- 
mia, para a Biologia, para a Quimica, para a Gestão de 
«stocks», etc, Para ilustrar a interpenetração do forma- 
lismo emergente da Comunicação noutras áreas, refi- 
ram-se, por exemplo, os estudos de qualidade de água, 
em que a mensagem a extrair é modulada (pela natu- 
reza) em parâmetros característicos como sejam, entre 
outros, o oxigénio dissolvido e a procura bioquímica de 


Teoria dos Fenômenos 


Electr. Din. Electr. Magn. Plasmas Acústica 


FIG. 4 


Teorias dos Fenómenos Fisicos Subjacentes 


oxigénio, e depois processada numa estação de trata- 
mento de águas. O homem como fonte, o homem como 
destinatário, desapareceu do ciclo de fluxo da informação, 

As figuras 4 e 5 ilustram alguns dos temas e disci- 
plinas abrangidos nas ramificações Teoria dos Fenóme- 
nos e Tecnologia consideradas na figura 2, Às listagens 
apresentadas privilegiam ramos de conhecimento tradi- 
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cionalmente ligados à engenharia electrotécnica. Com a 
diversificação das aplicações, a Teoria dos Fenómenos 
alarga-se para incluir, além de Ciências já referidas, ou- 
tras como a Meteorologia, a Economia, a Linguística, etc. 
Acentua-se que na figura 5 se consideram apenas tecno- 
logias referentes à implementação «hardware» dos siste- 
mas. A implementação «software» é aqui considerada 
como capitulo duma ciência nova, a Ciência dos Compu- 
tadores. 

Na figura 6 consideram-se disciplinas típicas no pro- 
cessamento, Como questão central apresenta-se a abs- 
tracção de modelos e suas representações formais. 


Sist.Mec. 


T.Quant. T.Quant Mat. Elect. Diel. 
dos do 
Solidos Magn 
(Metais, 
e ira 
Isolad.) 
FIG. 5 
Tecnologia 


No âmbito do formal, numerosos são os tópicos que, pelo 
seu desenvolvimento actual, se autonomizaram em disci- 
plinas, de elevado grau de abstracção e formalismo, como 
a Optimização, a Programação Matemática, a Identifica- 
ção, etc, Estas disciplinas concorrem em domínios fun- 
damentais das aplicações cuja nomenclatura é mais tra- 
dicional, como o Projecto de Circuitos, o Controlo, a 
Comunicação, etc. Estes incluem na figura 6, cami- 
nhando para a especialização e mais perto do real, dis- 


Computadores | +. (Programação) e-—— Algoritmos 


ciplinas tipicas como sejam no Controlo, os Servomeca- 
nismos, o Controlo Óptimo, a Investigação Operacional e 
na Comunicação, a Teoria Estatística da Comunicação, 
os Sistemas de Comunicação, etc. O suporte matemático 
encontra-se em ramos como a Análise (Variável Real, 
Variável Complexa, Funcional, etc.), a Teoria das Probabi- 
lidades e Processos Estocásticos, a Álgebra (Linear, Mo- 
derna, etc), as Matemáticas Aplicadas (Cálculo Opera- 
cional, Análise Numérica, Equações Diferenciais, etc.). 

Há tópicos que, mercê de circunstâncias várias, em- 
bora interpenetrando os desenvolvimentos anteriores, se 
especializaram no seu âmbito, nas suas preocupações, no 
seu formalismo, de tal modo que merecem uma referência 
autónoma divergente. Na figura 6 menciona-se, por exem- 
plo, o Processamento Digital, Esta disciplina, recente (um 
dos primeiros livros a ela dedicada data de 1968, [4], 
constitui hoje campo de actividade fértil que atingiu ele- 
vado grau de maturação revelado quer pelas revistas 
técnicas especializadas, p.e. [5], quer pelos numerosos 
livros já publicados, p.e. [6], [7], [8]. 

Um tópico fundamental no processamento é consti- 
tuido pelos algoritmos e pela sua teoria, Os algoritmos 
representam a solução do problema, não como uma fór- 
mula analítica (de significado nem sempre transparente), 
mas como uma rotina de operações que a exprimem em 
termos dum conjunto de números (tabela, gráfico), Os al- 
goritmos, uma vez programados num computador, tradu- 
zem uma implementação, da acção projectada, dum tipo 
novo — a implementação software. Do ponto de vista con- 
ceptual, e em termos do utilizador, esta implementação 
atinge os mesmos objectivos que podem ser atingidos 
pela implementação hardware, noutros contextos O utili- 
zador pode adoptar em relação a ela a mesma atitude, 
de caixa negra, que adopta muitas vezes em relação a 
uma implementação hardware, nomeadamente a atitude de 
se saber servir dela, saber utilizá-la, mesmo que não co- 
nheça, em pormenor, o seu funcionamento interno, À Teo- 
ria dos Algoritmos tem conhecido um incremento sempre 


[Processamento] 


Processamento Digital 


Tópicos Especializados 


Modelos Representação 
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Tópicos no Dominio do Processamento 


TÉCNICA 455 


171 


crescente, incluindo tópicos como a sua representação e 
implementação (fluxogramas, estruturas de controlo, sin- 
cronia e assincronia, etc.), algoritmos de natureza aritmé- 
tica (multiplicação, factorização com inteiros e polinó- 
mios, teste de números primos, avaliação de máximo divi- 
sor comum de polinómios), métodos matriciais, resolução 
de equações algébricas e diferenciais. Na escolha do al- 
goritmo a implementar há que atender a factores como a 
convergência, erros associados, memória ocupada, etc, 

A figura 6 ilustra a relação, através dos algoritmos, 
considerados como implementação software das acções 
programadas, entre o processamento e a Ciência que é, 
talvez, o simbolo dominante da sociedade tecnológica — a 
Ciência dos Computadores, Os computadores, que-vão 
buscar as suas raizes nos anos 40 a problemas bem con- 
cretos e a necessidades bem reais das aplicações, aban- 
donam a fronteira difusa do interface com as outras 
ciências e constituem-se num ramo próprio em expansão 
complexa e acelerada, ver a figura 7. Como imagem com- 


Ciência do Computador 
Software ao Hardware 
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Matemáticas: Lógica, Estruturas Algêébricas,Com- 
binatória,... 


FIG. 7 


Ciência do Computador 


pare-se a génese de novas ciências, a partir de uma 
ramificação comum, à multiplicação, por cissiparidade da 
célula, na Biologia, Este fenómeno aliás continua a ma- 
nifestar-se como é exemplo recente a engenharia de 
software, 

Na figura 7 desdobrou-se a Ciência do Computador 
em dois ramos básicos, No ramo de Hardware incluem-se, 
por exemplo, problemas como: i) estruturas de compu- 
tação, l.e., arquitecturas do computador para projecto de 
hardware capaz de suportar determinadas características 
de software desejadas; ii) sistemas de informação em que 
sobressaiem problemas relacionados com o objectivo de 
acesso a uma comunidade de utilizadores com as ine- 
rentes limitações técnicas (sistema de memória, sincroni- 
zação, fiabilidade, etc.); lil) estruturas de dados, técnicas 
gráficas, etc., etc. No ramo de Software incluem-se a pro- 
gramação (com a sua interacção e compromissos com o 
hardware), as linguagens, com a sua sintaxe e semântica, 
as técnicas de geração de códigos, a representação de 
sistemas de dados, etc, Evidentemente que esta divisão 
não é estanque, devendo explorar-se a relação entre a 
linguagem de programação, a estrutura dos dados, o pro- 
jecto da máquina por forma a optimizar a interacção 
entre a estrutura algoritmica, a linguagem, a estrutura de 
dados e o hardware. 
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A Ciência do Computador, e os seus desenvolvi- 
mentos especificos, como sejam a inteligência artificial, 
os sistemas de programação automática, as linguagens 
simbólicas, etc., fertilizaram campos da Matemática como 
a Álgebra (Sistemas Algébricos, Teoria dos Conjuntos, 
Homologia, etc.) a Teoria dos Números, a Combinatória 
(Optimização Discreta, a Teoria dos Grafos, Programa- 
ção Linear e Inteira, Jogos, etc.) a Lógica (Sistemas 
Formais), as Teorias da Computabilidade e Recursibili- 
dade, etc. 

A arrumação anterior não pressupõe que o conheci- 
mento seja estanque, emparedado por barreiras rigidas. 
ântes pelo contrário. Os diferentes ramos de conheci- 
mento entrelaçam-se numa teia em que por vezes os nós 
são de distinção dificil. A figura 8 exemplifica alguns sis- 
temas que ilustram a interface entre os ramos anteriores. 
Por exemplo a antena, enquanto radiador (ou captador) de 
energia interpenetra com o Electromagnetismo, com a Pro- 
pagação e a Radiação. Mas a antena, enquanto processa- 
dor de informação, é tópico privilegiado na Teoria da 


Electr.Magn. Antenas 


(a) 


Computadores Sist.Dig. 


| 


Process.Multidim. 


Electrônica 
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Frocessamento 


Semiconductores 


(tc) 
FIG. B 


Interfaces: Exemplos 


Comunicação. Nos sistemas digitais concorrem os compu- 
tadores, o processamento, a electrónica. Nos sistemas 
de comunicação óptica convergem teorias das ondas su- 
perficiais (Electromagnetismo, etc.), técnicas frequencialis 
de processamento, tecnologias de dispositivos semicon- 
dutores dependentes de processos de efeitos acoplados 
(piezoelectricidade, magnetoelasticidade, etc.). 

Neste relatório identifica-se a Teoria dos Sistemas 
(TS) com a sua ramificação processamento, tal como 
explicitada na figura 6. À TS é assim um divisor comum 
a várias aplicações. Plagiando a expressão (usada noutro 
contexto) dum professor dos cursos do ano propedêutico, 
a TS é daltónica — conserva dos problemas o que há 
de mais essencial (modelo), formaliza-os em termos ma- 
temáticos (representação) e analisa estas imagens do 
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real, projectando novos sistemas, novos sinais. Mas a TS 
só é fecunda se atender à natureza intrinseca dos fenó- 
menos subjacentes e quando a implementação dos seus 
projectos é testada na confrontação com o real. 


2. Três Questões Básicas na Teoria 
dos Sistemas 


Na secção anterior procurou-se estabelecer o âmbito 
da Teoria dos Sistemas (TS) referindo, através de uma 
listagem exemplificativa, de modo nenhum completa, do- 
mínios concretos por ela abrangidos. Da análise dos 
diagramas apresentados importa abstrair o ponto de vista, 
ou seja, a metodologia subjacente. Em poucas palavras 
a TS corresponde ao desenvolvimento formal de técnicas 
comuns a uma grande diversidade de problemas. 

Estas técnicas dizem respeito a três questões: 


|) Modelo/Representação do Sistema. 
ii) Análise do Sistema. 
jii) Síntese de Sistemas. 


Desenvolve-se brevemente cada um destes pontos. 


2.1 Modelo/Representação 


Qualquer que seja a aplicação de interesse prático, 
que designaremos genericamente por sistema, a TS esta- 
belece hipóteses de comportamento que, tomadas no seu 
conjunto, se referem por modelo do sistema. O modelo é 
a entidade, imagem do sistema, que interessa à TS. Dife- 
rentes hipóteses corresponderão a diferentes modelos. 
Reciprocamente o mesmo modelo serve de abstracção a 
diferentes sistemas. O modelo representa um compro- 
misso entre o real e o formal, compromisso este ditado 
por factores de natureza consciente (e por isso mesmo 
controláveis pela Teoria), decorrentes da necessidade de 
simplificar o sistema complexo para o sujeitar ao estudo, 
ou de natureza inconsciente, como resultado de ignorân- 
cia da fenomenologia real. 

A passagem do sistema (real) para o modelo (formal) 
corresponde hoje em dia, com a diversificação das aplica- 
ções, a um dos domínios mais importantes de TS. Pode 
designar-se genericamente por Modelação. A experiência 
adquirida em diversos campos levou à sistematização de 
modelos em classes distintas como por ex., modelos ll- 
neares ou não lineares, variáveis no tempo ou invarian- 
tes, continuos ou discretos, determinísticos ou probabi- 
lísticos, parâmetros concentrados ou distribuídos, etc., ver 
pe. [2]. 

Esta classificação dos modelos (ela também sempre 
incompleta, podendo ser mais fina ou mais grosseira) 
correspondeu (e deu origem) a desenvolvimentos mate- 
máticos específicos. À classe mais rica de resultados é 
sem dúvida a dos modelos lineares, invariantes no tempo. 

O segundo aspecto que importa comentar neste pa- 
rágrafo é o da representação dos modelos, ou seja, a 
tradução das hipóteses simplificativas que se assumem 
na linguagem formal da matemática. Esta descrição ma- 
temática — a representação ——- depende de factores vá- 
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rios, como a escolha das coordenadas, do sistema de 
eixos, etc., sendo possível diferentes representações para 
o mesmo modelo. 

De entre todas as possiveis representações reportar- 
-nos-emos a duas que mereceram grande desenvolvimento 
na TS e caracterizaram periodos distintos da sua evolu- 
ção. Referimo-nos à representação exterior (Rex) e à re- 
presentação interior (Rin).A discussão prossegue no con- 
texto de modelos lineares, invariantes no tempo. 

Rex, com raizes na Teoria dos Servomecanismos 
Clássicos (Hazen, Nyquist, Bode, Evans) e na análise 
frequencial do problema da Comunicação (Nyquist, Hart- 
ley, Black), descreve o sistema por uma função de trans- 
ferência e estuda o seu comportamento no domínio dual 
do tempo — o chamado espaço da frequência (com- 
plexa). O instrumento de análise é o Cálculo Operacio- 
nal (Transformadas de Laplace, Fourier e Z). Constitulu 
a representação privilegiada nos anos do pós-guerra, até 
final da década de 50. 

Rin, particularmente adaptada ao estudo dos proble- 
mas surgidos com os programas espaciais, caracteriza O 
sistema por um conjunto de equações diferenciais ordi- 
nárias de 1.º ordem, estudando-o directamente no domi- 
nio do tempo. O instrumento de análise é a Teoria das 
Equações Diferenciais. Dominou os avanços na TS nos 
anos 60. 

A 2.º representação permitiu melhor conhecimento 
da estrutura dos problemas básicos na TS. Esta repre- 
sentação, aliada ao (e motivado pelo) desenvolvimento 
e disseminação dos computadores (o primeiro grande 
computador transistorizado aparece em 1958) substitui ao 
ponto de vista analítico o ponto de vista algoritmico j.e., 
a solução do problema é encarada, não como uma fór- 
mula (quantas vezes sem significado aparente) mas sim 
como uma rotina de instruções facilmente programável 
e executável por um computador, apresentando a res- 
posta ao problema como uma tabela de valores numéricos 
ou um gráfico. Rin, à qual se associa o nome de Kalman, 
generaliza sem dificuldades à classe de sistemas multi- 
variáveis, os resultados anteriores da Teoria e, sobretudo, 
evidencia as limitações de Rex. O mérito de Rin residiu 
sobretudo em ter atraído a atenção para o estudo de 
pontos básicos essenciais como a controlabilidade, a 
observabilidade e outras questões relacionadas, o que 
permitiu grandes avanços teóricos. Nos anos 70, na se- 
quência de trabalhos, em que é justo citar o nome de 
Rosenbrock, retoma-se a representação exterior Rex, re- 
definindo em novos termos o operador representativo do 
sistema, O trabalho recente com Rex, retomando o Cál- 
culo Operacional generalizou aos sistemas multivariáveis, 
técnicas gráficas, de projecto, clássicas. Este regresso a 
Rex é no entanto um regresso algebrizado, i.e, é a 
introdução de conceitos algébricos na representação ex- 
terna. O instrumento analítico é agora o da Teoria das 
Matrizes Racionais (i.e. cujos elementos são funções 
racionais da variável complexa). 

No contexto actual pode dizer-se que para sistemas 
lineares invariantes no tempo, e dependente das aplica- 
ções ou da natureza dos problemas, há um ressurgimento 
de trabalho (teórico) no domínio da Rex. A Rex generaliza 


a uma classe mais vasta de problemas os resultados da 
Rin para esta classe de sistemas. 

Para sistemas lineares variáveis e sistemas não linea- 
res, Rin continua a revelar-se a via com maiores poten- 
cialidades. 

Notar que as técnicas de projecto, desenvolvidas 
quer no contexto de Rex quer no contexto de Rin, quando 
aplicadas aos problemas de interesse nas aplicações 
estão irremediavelmente ligadas aos computadores. 


22 Análise do Sistema 


A análise corresponde à necessidade de formaliza- 
ção da questão clássica das ciências físicas e da natu- 
reza, nomeadamente a de compreender o comportamento 
do sistema. Para evitar a descrição exaustiva da resposta 
a todas as possiveis (e diferentes) condições de funcio- 
namento, caracteriza-se o comportamento do sistema por 
um determinado número de propriedades e grandezas. 
No domínio dos sistemas lineares invariantes no tempo 
(SLIT), em que é válido o principio da sobreposição 
linear, assumem particular relevo grandezas como as fre- 
quências naturais e pólos, zeros, modos, etc. 


A análise procura pois conhecer as características 
essenciais das respostas sem ter que as calcular. Por 
exemplo, estabelecem-se as características de estabili- 
dade dos pontos de equilibrio, determina-se a possibili- 
dade, ou não, de controlar o sistema ou de reconstruir 
o seu estado (i.e. o conteúdo da sua memória) a partir 
de observações incompletas ou ruidosas, sem ter que 
fazer um estudo exaustivo de todas as situações. À aná- 
lise corresponde à determinação das propriedades estru- 
turais do sistema. 

Os trabalhos de Kalman sobre algebrização dos SLIT 
contribuiram para a compreensão da estrutura dos sis- 
temas. A questão essencial quase se pode resumir, em 
face de um objectivo determinado, a saber fazer a per- 
gunta certa. Estes trabalhos vieram a fertilizar muito 
significativamente os desenvolvimentos da TS ainda em 
outros dominios, com outros modelos, nomeadamente de 
parâmetros distribuidos, de dimensão infinita, lineares va- 
riáveis no tempo, certas classes de sistemas não linea- 
res (p. e. bilinear), etc. 

Um outro ponto que interessa realçar no dominio da 
análise e que caracteriza também a TS é a avaliação do 
comportamento do sistema por técnicas de simulação. 
O computador (digital) é a ferramenta base. O compor- 
tamento do sistema, na impossibilidade de ser testado 
na vida real (ou porque a sua natureza não permite a 
experimentação directa — p.e. sistemas sociais — ou por- 
que não é economicamente viável) é ensalado via simu- 
lação no computador digital. A simulação abriu campo 
extenso à aplicação destas técnicas e justifica a vanta- 
gem da caracterização do sistema por modelos mais com- 
plexos, como p.e. não lineares, 


2.3 Sintese 


A sintese ou projecto de sistemas (e sinais) é talvez 
o ponto básico que distingue a TS de outras ciências da 
natureza. Estas procuram (passivamente) compreender ou 
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explicar o comportamento do sistema enquanto, pela sin- 
tese, a TS intervém na sua evolução. 

Também na sintese o tipo de representação utilizada, 
Rex versus Rin, levou a que se desenvolvessem diferen- 
tes metodologias. Na sintese foi necessário ultrapassar 
a apreciação qualitativa do comportamento do sistema. 
à introdução de medidas quantitativas (funcionais de 
custo ou de qualidade) colocou o problema no dominio 
do Cálculo das Variações — é a Optimização. Assume 
particular relevo a classe dos SLIT com funcionais de 
custo quadráticas, uma vez que estas conduzem, em 
última análise, a um problema formal linear. 

Tipico de Rex são técnicas de implementação por 
factorização espectral e o processamento tipo batch 
(não recursivo), enquanto Rin veio a desenvolver técnicas 
algorítmicas de realização, caracteristicamente recursivas. 


3. O Departamento de Teoria 
dos Sistemas do IST 


Nas secções anteriores apresentaram-se uma radio- 
grafia (passe a expressão) (secção 2) e uma sistemati- 
zação das principais questões da TS (secção 3). Abor- 
da-se agora o problema do ensino da TS no curso da 
engenharia electrotécnica (EE). 

O ensino da TS enfrenta restrições típicas e anta- 
gônicas das quais se destacam: 


ij O formalismo e abstracção inerentes ao desen- 
volvimento das sofisticadas técnicas de modela- 
ção, representação, análise e sintese de siste- 
mas; 


li) À necessidade de introduzir os contextos bási- 
cos de aplicação mais característicos do dominio 
da EE, nomeadamente a produção, transporte e 
distribuição de energia, a comunicação, o con- 
trolo, a tecnologia, a instrumentação, etc.; 


iii) O ensino das metodologias de i) e dos contextos 
de il) ter que ser feito num periodo limitado 
(actualmente 10 semestres) com as necessárias 
extensão e profundidade; 


iv) À motivação dos discentes para que acreditem na 
teoria e simultaneamente desenvolvam o seu es- 
pirito critico em relação á(s) (insuficiências da) 
mesma teoria. 


Perante restrições contraditórias e dificilmente quan- 
tificáveis a solução tradicional tem sido a de apresentar 
as metodologias e técnicas formais no contexto das apli- 
cações a que os alunos são tipicamente mais sensíveis 
ou que sejam mais importantes do ponto de vista da sua 
vida profissional. Como exemplo, aponta-se a Teoria dos 
Circuitos que tem constituido nos curriculos das várias 
escolas de engenharia e ao longo do tempo e da evo- 
lução daqueles, o campo privilegiado de introdução des- 
tas técnicas formais. Por isso mesmo a Teoria dos Cir- 
cuitos é habitualmente a primeira disciplina de introdução 
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à conceptualização e ao ferramental analítico/algoritmico 
no contexto da EE. 

O Departamento de Teoria dos Sistemas (DTS) do 
curso de EE do IST faz incidir a sua actividade, quase 
se pode dizer exclusivamente, de acordo com a classifi- 
cação da figura 2 da secção 1, no ramo do processa- 
mento. Isto não significa que o DTS esgote os principais 
tópicos apresentados na figura 6 ou que no IST, pre- 
sentemente, os outros Departamentos não tenham também 
uma incidência maior ou menor nesta ramificação. De 
facto e até tendo em vista as restrições |) a iv) ante- 
riormente apontadas é de toda a vantagem que nos outros 
Departamentos, no contexto de aplicações especificas, se 
retomem, ou mesmo apresentem pela primeira vez, os 
métodos formais característicos da TS. 

O DTS inclui neste momento 13 disciplinas discrimi- 
nadas no quadro da figura 9, 

Para além das disciplinas desta figura interessa ter 
presente na discussão que se segue a disciplina de Fun- 
damentos de Telecomunicações (4º ano) do Departa- 
mento de Telecomunicações, onde se ministram aos alu- 
nos os conhecimentos básicos sobre Teoria da Informa- 
ção e análise (frequencial) de sistemas de comunicação. 
Outras disciplinas, situadas embora no âmbito do pro- 
cessamento, constituem a abordagem de tópicos inseri- 


Ano Disciplina Conteúdo 

Pa aist. Lógicos Alg. Boole, Máq. Combin. e 
sequence. 

ds Computadores Organização dum Compu- 
tador 

3.º |Mat. Apl. |, | Var. Compl., Cálc. Op., F. 
Gener, 

o Ra Circ. Comp. Introdução à Teoria 
dos Circuitos 

e Lab. Intr. Laboratório de Apresenta- 

 Circ. Sin. ção de Fenomenologia de 

interesse em EE, com ênfa- 
se na Comunicação 

3.º 1|T. Sist. Sin. Representação de Sin,, 
Lapl., Four., Convolução, 
Quadripolos 

4º Sist. Dig. |, II Implementação de Circ. 
com tecnologia digital 

4º | Servocom. |, I| | Controlo Clássico (SLIT, 
Unidim., F. Transferência) 

5. Servo. Ill Controlo Clássico 
(complementos, Projecto) 

5.º aist. Contr. Sist. Dinâmicos Multidim. 


— técnica de projecto em 
Contr. Ópt. e Comun. Est. 


Comun. 


FIG. 9 


Departamento Teoria dos Sistemas 
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dos em aplicações especializadas pelo que não serão 
explicitamente consideradas na sequência. 

Das 13 disciplinas do quadro, e comparando-as com 
o diagrama da figura 6, conclui-se que Sist. Lógicos, 
Computadores e Sistemas Digitais | e Il formam um 
núcleo (ainda elementar) especializado do DTS num do- 
minio de interface com a Ciência dos Computadores. 

As Matemáticas Aplicadas | e Il constituem discipli- 
nas de aprofundamento de ramos da Matemática com 
particular incidência na EE em geral e na TS em parti- 
cular, devendo abranger metodologias necessárias às di- 
ferentes ramificações desta. 

As restantes sete disciplinas (bem como algumas 
outras de que se recorda a já referenciada Fundamen- 
tos de Telecomunicações) constituem um conjunto onde 
se exploram, de pontos de vista muitas vezes comple- 
mentares, assuntos e preocupações afins. 

Estas disciplinas não esgotam os campos discrimi- 
nados, explicita ou implicitamente na figura 6. E escolhem 
para contexto especifico de apresentação dessas meto- 
dologias a Teoria dos Circuitos, o Controlo e a Comu- 
nicação. 

A disciplina de Circuitos e Computação é de natu- 
reza formal, constituindo actualmente divisor comum de 
todos os alunos da EE. É uma introdução, no contexto 
dos circuitos, a metodologias de análise (tanto do do- 
minio do tempo como no domínio da frequência) e a 
técnicas elementares de sintese. Pde em realce os pon- 
tos de vista analítico e algoritmico, procura sensibilizar 
os alunos simultaneamente para o rigor do formal (via 
axiomática) e para as limitações da Teoria, introduz a 
implementação hardware (com alguma prática laborato- 
rial) e software (simulação). Sendo a primeira aborda- 
gem dos alunos a estas matérias constitui, passe a ex- 
pressão, a primeira sintese necessária ao processo do 
conhecimento. 

Na Teoria dos Sistemas e Sinais aprofunda-se, no 
contexto dos Circuitos e das Comunicações, a interacção 
entre os sistemas e sinais, abordando-se a representa- 
ção frequencial de sistemas e sinais bem como técnicas 
mais especializadas de análise. 

Os Servomecanismos |, Il e |ll prosseguem presen- 
temente a chamada via clássica do controlo, concen- 
trando a sua acção em sistemas unidimensionais, linea- 
res, invariantes, desenvolvendo a representação por fun- 
ção de transferência, com incidência formal no Cálculo 
Operacional, privilegiando-se técnicas de projecto basea- 
das nos métodos de Nyquist, Bode, Evans. Os Servome- 
canismos Ill terão tendência a transformar-se numa dis- 
ciplina de projecto e instrumentação. 

A disciplina Laboratório de Introdução aos Circuitos 
e Sinais corresponde a uma conceptualização, por via 
heurística e experimental, de fenomenologia da EE, com 
particular incidência no domínio da Comunicação e do 
Controlo. Pretende constituir uma motivação ou ponto de 
referência para outras disciplinas de natureza formal. 
Defende-se a inserção de disciplinas com este cariz em 
pontos mais recuados do curso (2.º semestre do 1.º ano, 
1.º semestre do 2.º ano). Deve notar-se que a experi- 
mentação rigorosa nestes domínios, em particular no Con- 
trolo e na Comunicação, é para ser desenvolvida noutras 
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disciplinas, quer actualmente já existentes, quer possi- 
velmente a criar. Há necessidade de ultrapassar neste 
domínio as conhecidas dificuldades de ordem material 
(instalações, equipamento, etc.) e também humana. 

É evidente que o panorama que se tem vindo a 
descrever como abrangido pelas disciplinas mencionadas 
deixa vastos campos, no ramo do processamento, ainda 
por cobrir. A dificuldade de situar a disciplina de Sis- 
temas de Controlo e Comunicação não reside em encon- 
trar o que ensinar, antes pelo contrário, reside na escolha 
dos tópicos a abordar. 


4. Sistemas de Controlo e de Comunicação 


Discute-se em primeiro lugar a própria designação 
escolhida. Ela junta os três vocábulos Sistemas, Controlo 
e Comunicação como resposta a duas preocupações: 


ij) & da abstracção formal e universalidade dos con- 
ceitos abordados, inerentes à Teoria dos Siste- 
mas, tal como se tem vindo a explicitar. 


ii) À da exploração dos contextos particulares, (Con- 
trolo e Comunicação) hoje em dia ainda os mais 
significativos de aplicação da TS. 


Uma disciplina destas tem necessariamente que abor- 
dar o seu programa e o seu conteúdo estruturar-se-fo em 
torno dos três pontos já referidos: i) Modelo/Represen- 
tação; il) Análise iii) Sintese, Esta sistematização não 
significa uma abordagem sequencial destes três tópicos, 
que se revela má pedagogicamente. A necessidade de 
motivação, e o alertar para as próprias dificuldades ine- 
rentes ao real, aconselham a que se retomem ao longo 
do tempo os três pontos anteriores, abstraindo novos 
modelos ou aprofundando outras técnicas formais mais 
complexas. 


As outras disciplinas em que TS é abordada de 
forma sistemática (ver a discussão na secção 3), circuns- 
crevem o seu âmbito quase exclusivamente a siste- 
mas lineares, invariantes, deterministicos, unidimensionais. 
A disciplina de Sistemas de Controlo e de Comunicação 
(SCC) tem pois, à partida, 4 dimensões distintas em que 
pode progredir, consoante os sistemas estudados são: 


i) Lineares ou Não Lineares; 

ii) Invariantes (no Tempo) ou Variáveis (no Tempo): 
li) Deterministicos ou Probabilísticos; 

iv) Unidimensionais (') ou Multidimensionais. 
Pretende-se que SCC constitua uma visão de sintese 


nestes dominios, pelo que, em vez de se considerar ex- 
clusivamente o aprofundamento de uma ou outra destas 


hipóteses se optou por uma abordagem que considera 
as 4 dicotomias i) a iv). Esta opção é feita adoptando o 
quadro formal da representação interna, ver secção 3, 
o que permite que a mesma estrutura teórica estude 
simultaneamente as quatro direcções indicadas, sem ne- 
cessidade de reformular todo o contexto teórico quando 
se modifica uma ou mais das referidas dimensões num 
ou noutro sentido. O resto da secção considera as impli- 
cações na natureza de SCC dos quatro pontos i) a iv). 

Sobre a hipótese iv), com sistemas multidimensionais 
de dimensão finita, a dimensionalidade não restringe a 
natureza dos resultados obtidos, desde que o problema 
seja colocado no quadro teórico apropriado. A linguagem 
é a dos vectores e matrizes em espaços Euclideanos, 
sendo, para a maioria dos problemas no contexto de 
SLIT, a Álgebra Linear a estrutura formal conveniente. 
Como já se notou anteriormente, esta algebrização, moti- 
vada sobretudo pelos trabalhos de Kalman do final da 
década de 50, sugeriu também, no contexto da represen- 
tação externa, com os trabalhos de Rosenbrock, a gene- 
ralização dos métodos frequenciais clássicos ao problema 
dos sistemas multidimensionais. À partida pode pois afir- 
mar-se que SCC considera com toda a generalidade sis- 
temas multidimensionais. 

Na dimensão temporal, hipótese ii), conservando a 
linearidade, a generalização no sentido de sistemas, 
invariantes —» variáveis 
não faz normalmente perder a analiticidade do problema. 
Isto é, a teoria permite ainda, normalmente, exprimir os 
resultados em termos de uma expressão ou fórmula. 
Vejam-se como casos típicos a existência de uma matriz 
de transição ou os critérios de controlabilidade e obser- 
vabilidade. As consequências são de outra natureza, 
traduzindo-se na perda de aritmeticidade e de estaciona- 
ridade (ou uniformidade). Ou seja, as expressões anali- 
ticas a que se chegam, no caso de sistemas lineares, 
variáveis, são, habitualmente, complexas de interpretar, 
exigindo, a solução exacta, por via de regra, o recurso 
a algoritmos numéricos especializados com um número 
infinito de passos. Por outro lado, enquanto com sistemas 
invariantes as propriedades estruturais são independen- 
tes do instante em que se analisa o sistema, com siste- 
mas variáveis elas dependem, de um modo geral, do 
instante inicial. Conclui-se, do dito, que, ao conside- 
rar sistemas variáveis, a dificuldade do problema não 
reside na analiticidade dos resultados, mas sim na sua 
implementação e no teste dos critérios disponíveis. 
O esquema formal teórico é necessariamente mais sofis- 
ticado, com recurso mais frequente às subtilezas da 

análise, 

Sobre a linearidade versus não linearidade, dimensão 
|) supra, é de notar que os sistemas lineares se definem 
como a classe de sistemas que observa o princípio da 
sobreposição linear, enquanto que os sistemas não linea- 
res são genericamente caracterizados pela ausência dessa 
propriedade. Daí que o corpo de resultados disponiveis 


(1) Unidimensionais refere-se aqui a sistemas de entrada simples e saída simples. 
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seja no segundo caso muito pobre. Ao caminhar no 
sentido 

linear ——» não linear 
perde-se, por um lado, normalmente, a analiticidade e, 
por outro, ganha-se uma complexificação das questões, 
tendo que se considerar explicitamente o desdobrar do 
quadro teórico numa dimensão 


local versus global, 


Isto é, os resultados são referidos ao estado inicial (con- 
dição fronteira) e sua vizinhança. 

Finalmente, comenta-se a hipótese lii), dos sistemas 
determinísticos e dos sistemas probabilísticos. Os curri- 
cula dos cursos da EE, e em particular no IST, privile- 
giam de forma inequivoca a visão determinística do real. 
Há problemas em que, pela sua natureza intrinseca, esta 
visão limita significativamente a capacidade da teoria 
em compreender e influenciar correctamente a realidade. 
Nos sistemas probabilisticos assume particular importân- 
cia a classe dos processos Gaussianos e Markovianos. 
E isto porque, quanto aos processos Gaussianos: 


i) A linearidade preserva a Gaussianidade. 


ij) Os processos Gaussianos são completaments 
descritos por dois parâmetros estatísticos, a mé- 
dia e a covariância. 


ii) O Teorema do Limite Central justifica que nume- 
rosos processos ocorrendo na natureza possam 
ser considerados como assintoticamente Gaus- 
sianos. 


E, quanto aos procesos Markovianos: 


iv) Os processos Markovianos de ordem n são com- 
pletamente descritos em termos da sua densidade 
de primeira ordem e das densidades condiciona- 
das até à ordem n. 


Uma outra hipótese que simplifica a aritmética dos 
resultados é a da estacionaridade dos processos. 

Com respeito a sistemas probabilisticos e, tendo em 
conta as outras dimensões e a discussão anterior, pode 
dizer-se de forma sintética o seguinte: 


Não há acréscimo de 
dificuldades conceptuais 


Unidim. —> Multid. 


Fala-se em processos 
vectoriais 


Preserva-se Gaussianidade 


Inv. -—» Var. 


Perde-se Estacionaridade 


Perde-se Gaussianidade 
Preserva-se o carácter 


Markoviano 


Lin. => N. Lin. 


Perde-se Analiticidade 
Solução Algoritmica 
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Daqui resulta que, com SLIT, se associem habitualmente 
processos Gaussianos estacionários e, no outro extremo, 
com sistemas não lineares se assumam processos Mar- 
kovianos. 

Além das quatro dimensões discutidas poderiam con- 
siderar-se generalizações de outro tipo, como p.e., siste- 
mas de dimensão infinita. A principal questão refere-se 
aqui à natureza dos espaços envolvidos. Se a dimensão 
infinita for numerável há uma certa analogia com o que 
se disse sobre as dificuldades aritméticas na passagem 


Invariantes -> Variáveis. 


As questões fundamentais são essencialmente de natu- 
reza algoritmica. Se a dimensão infinita não é numerável 
a teoria reveste elevado grau de abstracção. 


Na disciplina SCC procura-se abranger no mesmo 
enquadramento teórico generalizações nas quatro dimen- 
sões i) a iv). Uma vez que se consideram sistemas mul- 
tidimensionais há necessidade de recorrer sistematica- 
mente a meios de cálculo automático no estudo dos sis- 
temas reais, como resposta natural aos encargos com- 
putacionais envolvidos. Toma-se pois a perspectiva de 
apresentar como solução prática das questões levanta- 
das a algoritmica, em prejuizo da solução analítica, Em 
contrapartida, a analiticidade é explorada no sentido de 
apresentar os conceitos no contexto mais simples, no- 
meadamente o de SLIT. Assim, o núcleo central do pro- 
grama é uma vez ainda constituido pelo estudo de SLIT, 
agora multidimensionais. Sem recorrer a uma necessária 
sofisticação dos argumentos da Análise Matemática, con- 
segue-se conter o formalismo dentro de limites razoáveis, 
apresentando os principais conceitos e métodos. Por 
outro lado referem-se, normalmente sem demonstração, 
os resultados correspondentes, para sistemas variáveis. 
Como exemplificativo da metodologia a adoptar no caso 
geral de sistemas não lineares, aborda-se um tópico nesse 
contexto. Escolhe-se a questão fundamental da Estabili- 
dade onde se apresenta com algum desenvolvimento a 
Teoria de Lyapunov. À razão da escolha reside no carác- 
ter formativo que se reconhece a essa teoria, não ape- 
nas em estudos de estabilidade como noutros campos. 


Antes da abordagem de técnicas de projecto em 
Comunicação introduz-se brevemente o ferramental ana- 
lítico necessário para a representação de processos esto- 
cásticos e para o estudo de sistemas prababilísticos. 

Em conclusão, a disciplina de Sistemas de Controlo 
e de Comunicação concentra-se em: 


1. Ji) Sistemas Lineares Invariantes Deterministicos 
Multidimensionais — o contexto matemático é 
o da Teoria das Equações Diferenciais Ordi- 
nárias de 1.º Ordem de Coeficientes Cons- 
tantes e o da Álgebra Linear. 


ij Os resultados teóricos são de natureza ana- 
lítica. 


li) Coloca-se a ênfase na realização em algo- 
ritmos numéricos, 
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iv) Incluem-se, por via de regra, sem demonstra- 
ção, resultados para sistemas lineares variá- 
veis. 


v) Os sistemas não lineares persistem em back- 
ground em SCC, Os sistemas lineares são 
apresentados como aproximação local dos sis- 
temas não lineares. 


2. Sistemas Probabilisticos — a ênfase é em siste- 
mas lineares, invariantes, multidimensionais e pro- 
cessos Gaussianos/Markovianos estacionários. 


3. Introduzem-se técnicas de projecto no Controlo & 


na Comunicação. 


5. Programa. Conteúdo. Métodos de Ensino. 
Bibliografia. 


5.1 Programa. Conteúdo. 
O esquema geral do programa de Sistemas de Con- 


trolo e de Comunicação divide-se, de acordo com a aná- 
lise feita anteriormente, em 


Parte | — Modelo/Representação 
Cap. | — Introdução à Teoria dos Sistemas. 
Cap. || — Fundamentos Matemáticos. 

Parte || — Análise 
Cap. Il — Sistemas Dinâmicos Lineares. 


Cap. IV — Controlabilidade e Observabilidade. 


Cap. V — Estabilidade. 
Cap. VI — Sistemas Discretos. 
Cap. VII — Sistemas Lineares e Processos Estocásticos. 


Parte Ill — Sintese 


Cap. VIII — Realização. 
Cap. IX — Compensação e Reconstrução. 
Cap. X 


Cap. Xl — Sistemas de Comunicação — Filtragem Li- 
near Óptima. 


— Elementos de Controlo Óptimo. 


Cap. XIl — Identificação. 


Este programa ignora numerosos tópicos, como p.e. 
a realização estocástica, o controlo sem interacção, o 
controlo estocástico, a inversão de sistemas, o controlo 


adaptativo, o «smoothing+, a predição, etc. Por outro 
lado, e como se verá, a abordagem dos campos enums- 
rados, designadamente na Parte Ill, é superficial, omitin- 
do-se pontos importantes. O esquema enunciado é um 
compromisso entre o tempo disponível e a vastidão do 
tema. Uma cobertura coerente exige necessariamente uma 
sequência de dois semestres. Refere-se seguidamente q 
conteudo de cada um dos capitulos. 


Nos dois capitulos da Parte | procura-se: |) motivar 
os alunos, apresentando-lhes uma perspectiva da evolu- 
ção da Teoria dos Sistemas (ver a Parte | de [2]) e dos 
principais tópicos a abordar; ii) explanar o ponto de vista 
adoptado; iii) justificar, na base de exemplos de aplica- 
ções concretas, o modelo e representação escolhidos; 
iv) apresentar os fundamentos matemáticos da Algebra 
Linear, da Teoria das Probabilidades e Processos Esto- 
cásticos e das técnicas de Optimização baseadas na Mi- 
nimização do Erro Quadrático Médio, necessários ao de- 
senvolvimento e compreensão da disciplina. 

A Parte || estuda as técnicas de análise e caracte- 
riza o sistema, introduzindo formalmente propriedades 
estruturais e os correspondentes critérios de teste. 

No cap. Ill considera-se a determinação analítica da 
resposta de sistemas lineares (variáveis e invariantes), 
a partir da matriz de transição, e as propriedades desta 
matriz. Desenvolve-se o cálculo da matriz exponencial 
(SLIT), por técnicas de Laplace e de Transformações 
de Semelhança. Segue-se o caso geral de sistemas não 
lineares, recordando-se técnicas de integração numérica 
de equações diferenciais ordinárias ('). O último tópico 
deste capitulo é a Linearização Dinâmica, ou seja, O 
método das perturbações para análise de sistemas não 
lineares, que permite gerar sistemas lineares variáveis a 
partir de sistemas não lineares. 


O capitulo IV formaliza em perguntas a fazer à estru- 
tura do sistema, questões como sejam: 


— é ou não possivel controlar o sistema? 


— é ou não possivel reconstruir o estado do sistema 
a partir das observações da saida? 


Introduzidos os conceitos de controlabilidade e de obser- 
vabilidade, desenvolvem-se critérios de teste para o caso 
simples de sistemas lineares invariantes. Os resultados 
referentes a sistemas lineares variáveis são também 
apresentados. O último tópico do capitulo refere-se à 
decomposição interna do sistema nas suas componentes 
controlável /observável, etc. 


O capitulo V desenvolve o importante tópico da Esta- 
bilidade. Após uma breve referência aos métodos clás- 
sicos, esclarece-se a natureza do problema, apresentan- 
do-se depois as definições formais de estabilidade de 
estados de equilibrio, no contexto mais geral de sistemas 
não lineares. Apresentam-se, para sistemas lineares, re- 
sultados baseados na natureza da resposta do sistema 
(1º método de Lyapunov). Finalmente considera-se q 


(1) No contexto de SCC trata-se apenas de referir a sua existência já que esses métodos são estudados na Análise Numérica e reto- 
mados, presentemente, na disciplina de Circuitos e Computação, 3.º ano de EE do IST. 
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Método Directo de Lyapunov (2.º método). Discutem-se 
brevemente algumas técnicas (algoritmicas) de geração 
de funções de Lyapunov para sistemas não lineares. 
Os resultados da teoria são particularizados ao contexto 
de sistemas lineares. 

O capitulo VI estabelece, para sistemas discretos 
no tempo, os resultados desenvolvidos nos capitulos II 
a V para sistemas continuos. 

O capitulo VIl, último da Parte Il, estuda a inter- 
acçlo entre sistemas e processos estocásticos. Admi- 
te-se, normalmente, que se conhece apenas a descrição 
incompleta de 2 ordem (limitação de natureza prática) 
dos processos. Dado que existe sempre um processo 
Gaussiano com a mesma descrição incompleta (i.e. com 
a mesma média e autocorrelação) consideram-se sobre- 
tudo processos Gaussianos, estabelecendo-se as equa- 
ções de propagação do 1. e 2. momentos. 

A Parte Ill estuda técnicas de sintese. No capi- 
tulo VIll é a realização de sistemas lineares invariantes, 
isto é, a obtenção da representação interna Rin (matrizes 
A, B, C, D) a partir da representação externa Rex (matriz 
de transferência), ou o problema inverso. As dificuldades 
residem em 


Po. ———e———————» Rin. 


Resolve-se o problema para sistemas com um nu- 
mero de zeros inferior ao número de pólos, realizando 
as chamadas representações canónicas. O caso de sis- 
temas de entrada simples/saida simples é estudado com 
maior desenvolvimento. O capitulo aborda ainda o pro- 
blema da transformação de coordenadas que permite, em 
certas condições, passar de uma representação genérica 
para uma das representações canônicas. 

No capitulo IX consideram-se simultaneamente pro- 
blemas típicos do controlo (compensação) e da comuni- 
cação (reconstrução), num contexto determinístico. No 
primeiro caso é o problema da colocação de pólos por 
realimentação linear negativa do estado e da saida. Este 
problema clássico do controlo é resolvido no contexto 
mais geral de sistemas multidimensionais. Na reconstru- 
ção considera-se o estabelecimento de condições que 
permitam a estabilização do observador. 


Os capitulos X e XI introduzem, respectivamente, 
no contexto do Controlo e da Comunicação, o projecto 
óptimo de sistemas. No capitulo X, enunciado o pro- 
blema, restringem-se as possíveis soluções pela escolha 
de uma funcional de custo que quantifica, no problema 
em estudo, as caracteristicas desejadas para o compor- 
tamento do sistema. Enuncia-se o Principio do Minimo 
de Pontryagin e apresenta-se a Programação Dinâmica 
de Bellman (no contexto discreto). Finalmente aplicam-se 
estes métodos ao projecto de sistemas lineares óptimos 
com diversas funcionais de custo (tempo minimo, com- 
bustivel mínimo, energia minima). 

No capitulo XI aborda-se o problema da estimação 
estatistica de formas de onda estocásticas. Descrito o 
problema no contexto dum sistema de comunicação (mo- 
dulação linear), classificam-se os estimadores e vários 
critérios de estimação. Adopta-se uma perspectiva Bayea- 
siana e apresentam-se os filtros recursivos lineares ópti- 
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mos de Kalmam-Bucy. Discutem-se ainda problemas de 
implementação. 

O capitulo XII, o último, é o da Identificação, cons- 
tituindo uma introdução a técnicas, alternativas às anali- 
ticas, de extracção de modelos a partir dos dados expe- 
rimentais, considerando-se métodos determinísticos e es- 
tocásticos. 


Note-se que, no esquema acima, não se reservou 
nenhum capítulo à apresentação exclusiva de algoritmos. 
Optou-se pelo desenvolvimento de alguns algoritmos no 
âmbito de problemas especificos nos diferentes capitulos. 


A discriminação completa dos capitulos e suas sec- 
ções é apresentada seguidamente. 


Parte | — MODELO /REPRESENTAÇÃO 


Cap. | — Introdução à Teoria dos Sistemas 


1. Perspectiva Histórica. 

2. aistema. Modelo. Representação. 
3. Exemplos de Aplicação. 

4. Sintese dos Resultados. 


Cap. || — Fundamentos Matemáticos 


1. Algebra Linear. 
Espaços Lineares. Transformações Lineares. Vec- 
tores Próprios. Valores Próprios. Funções de Ma- 
trizes. Teorema de Cayley-Hamilton. 


2. Elementos de Teoria das Probabilidades e Proces- 
sos Estocásticos. 
Espaços de Probabilidades. Independência. Proba- 
bilidades Condicionais. Lei de Bayes. 
Variáveis Aleatórias. Função Distribuição e Fun- 
ção Densidade de Probabilidade. Funções de Va- 
riáveis Aleatórias. Momentos. Média. Variância, 
Desvio Quadrático Médio. Exemplos. 
Variáveis Aleatórias Conjuntamente Distribuídas. 
Correlação. Covariância. Vectores Aleatórios. 
Função Distribuição. Função Densidade. Função 
Característica. 
Teorema do Limite Central. 
Processos Aleatórios. Caracterização Completa e 
de 2.º Ordem. Processos Aleatórios Conjunta- 
mente Distribuídos. Processos Aleatórios Vecto- 
riais. Função Autocorrelação. Função Correlação. 
Covariância e Covariância Cruzada. 
Estacionaridade. Ergodicidade. 
Função Densidade Espectral de Potência e Função 
Densidade Espectral de Potência Cruzada. 
Exemplos. Ruido Branco. Processo Gaussiano. 
Processo Markov. Processo de Incrementos Inde- 
pendentes. Processo Browniano. Processo Pois- 
son. 


3. Técnicas de Minimização do Erro Quadrático. 
Minimização do Erro Quadrático. Minimização Pon- 
derada do Erro Quadrático. Minimização Recur- 
siva, 
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Parte || — ANALISE 
Cap. Ill — Sistemas Dinâmicos Lineares 


1. Representação Interna, 
2. Sistemas Lineares Variáveis. Matriz de Transição. 
Série de Peano-Baker. Propriedades. 

. Fórmulas de Variação das Constantes, 

. Sistemas Lineares Invariantes. Matriz Exponencial, 
Propriedades. 

5. Cálculo da Matriz Exponencial. Transformada de 
Laplace. Transformações de Semelhança e Redu- 
ção a Formas Canónicas. 

6. Linearização Dinâmica. 

7. Algoritmos de Integração. 


+ Cu 


Cap. |V — Controlabilidade e Observabilidade 


1. Definição. 

2. Critérios para Sistemas Lineares Invariantes. 
3. Critérios para Sistemas Lineares Variáveis. 
4. Decomposição Interna. 


Cap. V — Estabilidade 


1. Introdução. 

2. Breve Referência aos Métodos Clássicos. 

3. Questões Preliminares. 

4, Estabilidade de Estados de Equilibrio: Definições. 
5. Estabilidade de Sistemas Lineares. 

Sistemas Lineares Variáveis no Tempo. 
Sistemas Lineares Invariantes. 

segundo Método ou Método Directo de Lyapunov. 
Introdução. Sistemas Não Lineares Invariantes. 
Sistemas Não Lineares. Geração de Funções de 
Lyapunov. Sistemas Lineares. 


Edo 


Cap. VI — Sistemas Discretos 


1. Sistemas Lineares Variáveis. Matriz de Transição. 

2. Sistemas Lineares Invariantes. Cálculo da Matriz 
de Transição. Transformada Z e Transformações 
de Semelhança. 

3. Fórmula de Variação das Constantes. 

4. Controlabilidade e Observabilidade. Critérios. 

5. Estabilidade. 


Cap. VIl — Sistemas Lineares e Processos Estocásticos 


1. Resposta do Sistema a Entradas Estocásticas. 

2. Equações de Propagação da Média e Covariância. 

3. Sistemas Lineares Invariantes. Regime Estacio- 
nário. 


Parte Ill — SINTESE 
Cap. VIII — Realização 


1. Matriz de Transferência e Representação Interna. 
2. Representações Canónicas. 
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3. Transformações de Coordenadas e Representa- 
ções Canónicas. 

4. Cancelamento de Pólos e Zeros. Efeito na Con- 
trolabilidade e Observabilidade. 


Cap. IX — Compensação e Reconstrução 


1, Colocação de Pólos. Realimentação Linear Nega- 
tiva do Estado. 
Sistemas de Entrada Simples. Sistemas de En- 
trada Múltipla. 

2. Reconstrução do Estado. Observadores. 

3. Colocação de Pólos. Realimentação Linear Nega- 
tiva da Saida. 

4. Observadores de Ordem Minima. 


Cap. X — Elementos de Controlo Óptimo 


1. Enunciado do Problema. Funcional de Custo. Res- 
trições no Controlo. Exemplos. Tempo Mínimo. 
Combustivel Minimo. Energia Mínima, 

2. Principio do Minimo de Pontryagin. 

3. Programação Dinâmica de Bellman. 

4. Projecto de Sistemas Lineares Óptimos com Custo 
Quadrático. 


Cap. XI — Sistemas de Comunicação — Filtragem Linear 
Óptima 


1. Enunciado do Problema. Estimadores (não) Pola- 
rizados, de Variância Mínima, Consistentes. 

2. Filtragem Recursiva Linear Óptima. 
Filtros de Kalman-Bucy. Equação de Propagação 
do Filtro. Equação de Ricatti da Covariância. 

3. Implementação. Algoritmos e Considerações de 
Natureza Numérica, 


Cap. XII — Identificação 


1. Métodos Determinísticos. 

2. Métodos Estocásticos. 
Técnicas de Autocorrelação. Análise por Séries no 
Tempo. 


Esta lista de assuntos, constituindo o núcleo conside- 
rado minimo indispensável numa disciplina com os objec- 
tivos de SCC, revela-se demasiado extensa para ser co- 
berta integralmente num único semestre. Discutem-se por 
isso programas restritos (alternativos) a prosseguir num 
semestre. 

Estabelecem-se, no entanto, primeiramente, dois pon- 
tos a ter em consideração: 


ij O semestre é dimensionado de acordo com o mé- 
todo exposto em [9], o que se traduz por 14 semanas 
Úteis de aulas. É de notar que, embora no passado ainda 
se não tenha atingido este total, é de esperar que, futu- 
ramente, isso possa suceder, com a adopção, no IST, de 
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calendários escolares e métodos pedagógicos mais ra- 
cionais. De qualquer forma programa-se, na sequência, 
a disciplina para 


14x2=28 sessões teóricas 
14x 1=14 sessões práticas 


em que cada sessão corresponde a um módulo de 1 2 h. 


ii) Os alunos que frequentaram até ao presente (1.º 
semestre 1977/78 inclusivé) as sucessivas versões de 
SCC evidenciaram grandes dificuldades com a Algebra 
Linear (para não falar da Teoria das Probabilidades e 
Processos Estocásticos). No entanto, os futuros alunos 
terão já frequentado a disciplina de Circuitos e Compu- 
tação (3.º ano de EE no IST), onde se faz extensivo uso 
da Teoria das Matrizes e se desenvolve, no contexto 
de circuitos, a representação interior. Espera-se que algu- 
mas das dificuldades até agora encontradas ,decorrentes 
da falta de destreza no cálculo matrical e da desmotiva- 
ção (por dificuldade em idealizar sistemas físicos — p.e. 
circuitos eléctricos de simulação analógica), venham a 
estar naturalmente resolvidas à partida. Pensa-se assim 
que os futuros alunos já terão tido contacto anterior com 
o ponto de vista adoptado em SCC. 


1 (Repres./Modelo) (3) 
1-1 (Alg. Linear) (0) 
III (Fesposta em Sist. Lin.) (6) 
V (Estabilidade) (4) 
Iv A4Control./Observ.) (2) 
| VI (Sist. Niscretos) (2) 
vem (Realização) (1) 
E, (Compensação e Reconstr.) (3). 
ati (tTêcn. de Min. de Err. Quad.) (1/2) 


NE II=2 (Pr. e Pr. Est.) (3) 
(6 1/2) VIT (Sist. e Pr, Est.) (1/2) 
(Contr.ôÓpt.) I 
XI (Sist. Comun. Filtr.- (3) 
A Lin.) 
XIT (Tdentificação) - po (0) 
TOTAL (28) 
FIG. 10 
DE [ad E 


Programações Exequíveis para SCC 


Tendo em atenção estés pressupostos propõem-se 
as programações (alternativas) da figura 10: Em parên- 


tesis indica-se.o número de sessões dedicado a cada 
tópico. Nota-se que a sequência apresentada anterior- 
mente foi alterada. Por outro lado as duas versões da 
figura 10 optam entre o Cap. X e o Cap. XI. Note-se 
que esta última inclui também os elementos de descrição 
probabilistica. Devido talvez às inclinações do autor deste 
relatório, e à sua preparação anterior, tem-se preferido 
esta versão à versão que opta pelo Cap. X. No futuro 
pensa-se alternar entre uma e outra via. É de notar que 
ao Cap. XII não se reserva explicitamente nenhuma sessão 
teórica, devido à falta de tempo. Breves referências ao 
problema poderão ser feitas no contexto do Cap. XI. 
Também os resultados mais sofisticados da Algebra Li- 
near e da Teoria das Matrizes (Cap. 1l-l) são introduzidos, 
quando necessários, no decorrer do Cap. Ill, enquanto os 
que se consideram mais elementares são recordados nas 
sessões de problemas. 


5.2 Métodos de Ensino 


Do que se disse no parágrafo anterior se conclul 
que a disciplina SCC é de natureza formal, adoptando-se 
a distribuição 2 do quadro 4.1 de [9] e a filosofia expla- 
nada em 3.2 de [10]. 


O seu ensino é assim feito a partir de: 


i) Sessões teóricas (duas por semana) de 1 % hora 
cada, onde se apresentam os pontos básicos da matéria; 


li) Sessões de problemas (uma para duas sessões 
teóricas) de 1/4 hora cada, onde se exemplificam com 
problemas, tanto quanto possível reflexo de aplicações 
de interesse, os conceitos e resultados da teoria ('). 


ii) Séries de problemas, de escopo limitado, a re- 
solver pelos próprios alunos, constituindo oportunidade, 
para estes, de estudo e reflexão sobre a matéria. 


Importa sobre este tema fazer dois comentários adi- 
clonais: 


) O primeiro refere-se à não consideração expli- 
cita, presentemente, no contexto da disciplina, de téc- 
nicas de simulação digital e analógica. Este ponto poderá 
vir a ser colmatado no futuro. De facto, e no contexto 
de Circuitos e Computação, está a desenvolver-se um 
«pacote» de subrotinas de aplicação generalizada, a que 
os alunos de SCC poderão também vir a ter acesso, 
Quanto à simulação analógica, o problema está relacio- 
nado com um melhor apetrechamento dos Laboratórios 
de Circuitos e Computação e de Servomecanismos. Tal 
como já se apontou anteriormente, as dificuldades não 
são apenas de natureza material. 
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li) O-segundo diz respeito à inserção de SCC no 
curriculum de EE (licenciatura). O seu presente isola- 
mento no contexto do curso poderá ser ultrapassado por 
duas vias. À primeira; inserindo já -estes-assuntos em 


ha so “holl = À. ê atri 5” 5) * fe 


(7) Nas versões anteriormente praticadas, e uma vez que o ensino tem sido exclusivamente ministrado pelo autor, as sessões são 
mais de natureza teórico-prática. No futuro preferir-se-á o presente formato. * BRST: + Ea k 
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disciplinas do 3.º ou 4.º ano (!) e desdobrando no 5.º ano 
a SCC numa sequência de duas — uma no 1º e outra 
no 2.º semestre, Este desdobramento permitirá uma cober- 
tura mais extensa dos tópicos referidos na Parte III, bem 
como de outros aí omitidos. À segundo via diz respeito 
à criação de cursos de pós-graduação, incluindo cober- 
tura mais aprofundada destes campos. 


5.3 Bibliografia 


Para o programa desenvolvido anteriormente não 
existe, do conhecimento do autor, um livro que se possa 
adoptar como referência única. A dificuldade reside no 
facto de a literatura existente se dirigir a estudantes com 
preparação matemática e motivação para estes temas, 
muito diversas das que presentemente revelam os alu- 
nos do IST. Por outro lado, a literatura disponivel pres- 
supõe que estes assuntos sejam retomados em diferentes 
momentos do curriculum académico, com diferentes 
graus de sofisticação, e abrangendo diferentes tópicos, 
Por isso os livros que se encontram são, por via de re- 
gra, de âmbito muito restrito e especializado em relação 
ao programa descrito, pecando (por excesso ou por de- 
feito) ma dificuldade formal com que os assuntos são 
apresentados. 

A referência [11] é manifestamente introdutória, co- 
brindo apenas tópicos da Parte Il e do cap. Mill, en- 
quanto [12] dedica cerca de 50% do seu conteúdo à 
Algebra Linear. Dai que a cobertura que faz nos restantes 
capitulos sobre temas da Parte Il e Ill acima considerados 
seja superficial, remetendo para problemas o desenvolvi- 
mento de questões fundamentais. Do ponto de vista de 
extensão de tópicos abordados, grau de sofisticação (mo- 
derado), preço, [12] é sem dúvida uma referência acon- 
selhável. Sobre o Cap. XI, e de entre a vasta literatura 
existente, destaca-se em particular [13] pela sua acesse- 
bilidade e pela ênfase em aspectos de implementação. 
Neste livro encontra-se ainda um breve resumo (cerca 
de duas dezenas de páginas) dos elementos necessários 
da Teoria das Probabilidades e Processos Estocásticos. 
Sobre o Cap. X refira-se o tratado clássico [14], presen- 
temente algo desactualizado. Refiram-se ainda [15] sobre 
os Cap. IX, X-4, XI [16] e [17] sobre o Cap. V [18] e 
[19] sobre a Parte Il. Como resposta às dificuldades 
apontadas tem-se optado por ir coligindo apontamentos 
próprios sobre os assuntos abordados, ver [20]. 


6. Conclusão 


Neste relatório procurou-se primeiramente delimitar 
o âmbito da Teoria dos Sistemas, especificando alguns 
dos seus campos privilegiados de acção (secção 1) e 
discutindo algumas questõeôs básicas (secção 2). Reflec- 
tiu-se depois, na secção 3, sobre o Departamento de 
Teoria dos Sistemas do curso de Engenharia Electro- 
técnica no IST, departamento cuja actividade se circuns- 


creve ao dominio da Teoria dos Sistemas e em que se 
insere a disciplina de Sistemas do Controlo e Comuni- 
cação. Na secção 4 apresentou-se o ponto de vista 
tomado nesta disciplina, definindo-se a sua natureza glo- 
bal. Finalmente na secção 5 estudou-se a organização da 
disciplina, apresentando o programa, o conteúdo, os 
métodos de ensino praticados e comentando (brevemente) 
(alguma) (dja bibliografia. 
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In the early stages of the planning of a new exploi- 
tation, the error in ore reserve estimation is an important 
factor that must not be disregarded. As this error can be 
computed by geostatistical methods, the risk of the ven- 
ture can be calculated for different recognition patterns 
in a certain asumption of economical conditions. 

A case story is described for an open pit phosphate 
deposit in which grade/tonnage curves and their confi- 
dence intervals are related to an economic function for 
three drilling patterns and decisions are taken on the 
basis of cost and value of supplementary information. 


INTRODUCTION 


A phosphate deposit located in the State of Goias, 
Brasil, was recognized by a 100 x 100m regular grid of 
134 drillings. Pieces of core from these drillings were 
analysed for P;Os and CaO and grades in P;Os from 
apatite were computed for each im length sample. Blocks 
50 x 50 x 10m were built over all the deposit for selec- 
tion purposes in this stage of planning. For each of these 
units, a kriging procedure as described by Matheron 
(1970) was adapted in order to estimate the average 
grade and kriging variance. 

As the deposit is not yet completely recognized in 
depth, global ore reserve estimation is unneeded, but 
some decisions must be taken about the economics of 
the venture in the zone near the surface. 

For this zone, an open pit and a pilot plant were 
designed and strategies of exploitation must be com- 
pared. 


STRUCTURAL ANALYSIS 


The Variogram is the basic tool of Geostatistics. 
The mathematical expression of this function is: 


= —+ 
2x(h) = ([Y(u+h) — YWP) (1) 


+ 
Y(h) is the variogram, 


() Research Senior Assistants, CVRMUL — IST — Ave. Rovisco 


Nos estágios iniciais do planeamento de uma nova 
exploração mineira, o erro cometido na estimação das 
reservas é um factor importante que não deve ser des- 
prezado. O cálculo deste erro por métodos Geoestatisti- 
cos permite a execução de uma análise de risco fixando 
certas condições económicas. 

É descrito no artigo um exemplo concreto de apli- 
cação destas técnicas. Trata-se de um depósito de fos- 
fato explorado a céu aberto para o qual foram obtidas 
curvas tonelagem/teor com os respectivos intervalos de 
confiança. Estabeleceu-se uma relação entre estas cur- 
vas e uma função económica para três malhas de reco- 
nhecimento, tendo as subsequentes decisões sido toma- 
das com base num balanço entre o custo e o valor da 
informação suplementar. 


Y(u+h) and Y(H) are the values taken by the Regionali- 
> 
zed Variable in two points linked by the vector h. 

This function, as Matheron (1970) and David (1975) 
showed, depends on the sample volume where Y(L) is 
calculated, 

In the case of this deposit, the basic information 
includes bm spaced channel samples which must be used 
in computation of variograms for small scale study. À nu- 
merical deconvolution as described by David (1975) was 
performed to recover punctual variograms from channel 
samples and core drilling samples. 

An example of such a variogram in EW direction is 
shown in Fig. 1 and a theoretical model of spherical type 
was fitted. 


The basic parameters shown in Fig. 1 are: 


Co — nugget effect 
C —sill 
a — range 


The meaning of these parameters are discussed by 
Matheron (1970), David (1975), Journel (1975). 

In particular, the «range» is the distance for which 
the autocorrelation function vanishes. 


Pais, Lisboa, Portugal. 


Paper presented at 15th APCOM Symposium, Brisbane, Australia, 1977. 
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EW Variogram of accumulation in P;Os 


Variograms of the same type were computed for NS 
and vertical directions and a geometric anisotropy was 
detected (ranges of 345m and 37.5 m). 

The Regionaiizeu Variable used is the accumulation 
in P;Os (grade x thickness) from apatite. 


KRIGING 


For each block, the kriging technique provides an 
optimal estimation of the average grade and the variance 
of this estimation. Assuming a normal distribution for the 


average grade estimate Z, the real value Z of the grade 
in a block lies in the interval 
at the 95% probability level, 


2-2 L2C2+2%, (2) 


If one takes a set of n blocks, the confidence interval 
o* the average grade is: 


+ 20% 1 2% 
um ZE 
m m 


(3) 


where Z is the kriging estimate of the average grade of 
n blocks and cp is the average kriging variance of n 
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blocks. In other words, there is a risk of 2.5% that the 

average grade of a certain set of n blocks is lower 
0 2x 

than Z — 


Using all the available information, a block by 
block estimation for all units of the deposit was per- 
formed using a three dimensionnal kriging program that 


provides estimates Z and variances o. 

For details in kriging theoretical basis, see Matheron 
(1970) and for practical applications, see Journel (1975) 
and David (1975). 


GRADE/TONNAGE CURVES 


A small area in the deposit was chosen to build 
grade/tonnage curves and study the sensitivity of an 
economic function to the confidence intervals of those 
curves. 

The selected area is shown in Fig. 2. 

For each of the 445 blocks of this area, a grade/ 
[tonnage curve was computed on the basis of kriged 
values for grades. 

This curve can be used for purposes of planning 
at this stage as David (1972) suggested. 

When more information is available, namely, smaller 
exploitation units, loader machines to be used, blast holes 
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e Given Drilling Pattern 


FIG. 2 | 
Location of the blocks and drillings in the selected area 


geometry, etc., a short term planning must be done using 
the method developed by Matheron (1975) and applied 
by Marechal (1975). 

The grade/tonnage curve for this area is shown in 
Fig. 3 and the grade/average grade curve for the same 
area is shown in Fig. 4. 

lt is assumed that there is no error in the density 
of the ore and so the number of blocks with grade 
higher than a certain cutoff is proportional to its tonnage 
and no confidence intervals are assigned to this curve. 

On the other hand, confidence intervals for the ave- 
rage grade can be computed from 3 and the curve of 
Fig. 4 shows these intervals for the 95%, probability 
level. 


SUPPLEMENTARY DRILLING PATTERNS 


In the same area, the effect of a first supplementary 
campaign of seven drillings of about 30m length (the 
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FIG. 3 
Grade/tonnage curve 
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Grade/average grade curve with 959%, confidence intervals 


same as the given drillings) and a second campaign of 
ten more drillings was simulated. For location of these, 
see Fig. 5. 


- Jo Joel! 
O A 


e Given Drilling Pattern 
o First Supplementary Drilling Pattern 
O Second Supplementary Drilling Pattern 


FIG. 5 


Simulated supplementary drilling campaigns 


As kriging depends only on the variogram and geo- 
metry of samples used for estimation of each block, the 
variance of kriging for all the blocks in the area can 
be computed in both cases (first and second supplemen- 
tary drilling campaigns). 

in the zone of the interesting cut-off grades (2 to 
6%) the grade/average grade confidence intervals were 
computed and are shown in Fig. 6. 

If it is assumed that, for a certain cutoff, the mana- 
gement wants to take a risk of 2.5% of producing ore 


with an average grade lower than E 20x /m (see 3), 
the useful curves for planning the exploitation are 0, 1, 2 
in Fig. 6. 

If the first drilling campaign is performed, and for 
the same risk of 2.5 %, the minimum average grade fore- 


seen is higher than if no supplementary drillings is made 
(curves 1 and O in Fig. 6). 


RISK ANALYSIS 


Assuming a simple linear model for an economical 
function relating tonnage and average grade, the global 
Profit P for all blocks mined can be written as — 


P = T(gnv—0) (4) 
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FIG. 6 
Grade/average grade curves and confidence intervals for the 


given drilling pattern (0), first supplementary drilling cam- 
paign (1) and second supplementary drilling campaign (2) 
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where: 


T — is the tonnage of blocks mines (tonnes) 

Gn — is the average grade of blocks mined (%) 

v — is the selling price of 1 ton of product corrected 
for recovery in the ore dressing plant 

C — is the cost of mining of 1 ton of ore 


T is related to cut-off grade by the curve in Fig. 3 
and gm is related to the same variable by the curve in 
Fig. 4. 

Assuming a risk of 25%, the minimum Profit ex- 
pected for a certain cut-off can be calculated from 4 
using the lower confidence interval (curve O in Fig. 6) 
for computing gm — 20x /m. 


So, there is 2.5 % of probability of obtaining a Profit 
lower than Py: 


Pu = T[(9m — 20x/m)v — C] (5) 


If the first supplementary campaign is performed, a 
cost S, of drilling must be compared to the gain of pre- 
cision in estimation: 


Pu, = Tgm — 20x, [mv — C]—S, (6) 
the second drilling campaign leads to: 
Pu, = T[(gm — 20x,/m)v — C]— S; (7) 


where S, is the cost of drilling. 


AP, is the increment of the minimum (at a 2.5 % 
risk) profit expected by performing the first campaign: 


20x 20x; 
AP, = — 1 ih (9) 
« M m 


if áPy, is positive, it is interesting, in the assumptions 
made, to perform the supplementary drilling. 


AP, can be written as Pu, — Pu or 


20%| 20x) 
APu, =T ( e v— 8 (10) 


q: m 


and measures the interest of the second drilling cam- 
paign, compared with the first one. 


Values of APyu, and dPy, were computed for the 
area of Fig. 5, for the range of cut-off grades shown in 


Fig. 6. 
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Results are show in Table 1º and 22, 
TABLE 1 


Increments of minimum profit expected for a 2.5 % risk 
by performing the first drilling campaign 


Cut-off o o 
o Tonnage 20x 20x | APy, 
grade a 
2 21.665 2.18 36.73 
3 20.655 2.33 37.63 
4 19.241 2.85 44.34 
5 17.091 3.04 41.24 
6 13.787 3.15 32.93 
TABLE 2 


Increments from the first to the second drilling campaign 
expected in the minimum profit for a 2.5 % risk 


Second and third columns of Tables 1 and 2 were 
built from Fig. 3 and Fig. 6. 


CONCLUSIONS 


In the early stages of mine planning, the cost of 
supplementary information can be compared to the gain 
in precision of the estimations if it is assumed an allowed 
risk for the venture. 

In this case story it is necessary to conjugate the 
exploitation strategy (which depends mainly on the cut- 
-off grade) with the plant strategy (which depends mainly 
on the average grade). For a certain cut-off grade, it can 
be foreseen the expected average grade from the curve 
ot Fig. 4. 


(2) For proprietv reasons values are multiplied by a constant. 
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The Profit is a function of this average grade, the 
lower limit of which may be fixed by a maximum allowed 
risk. 

Supplementary drilling provides narrower confidence 
intervals for average grade curves and so, for a fixed 
risk, the minimum profit expected will be higher. 

Drilling patterns used in this case story are arbitrary 
and so conclusions are empirical but a step-by-step pro- 
cedure may be attempted to optimize the method. 
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NOTICIÁRIO 


(Continuação da pág. 167) 


NO FINAL DO SÉCULO O SOL SERÁ CEM VEZES MAIS BARATO 


Afinal a crise da energia que tanto apregoam aos nossos ouvidos, deve ser bastante suave e limitada, pois 
que nenhum Governo aplicou, ainda, uma política de energia séria e incisiva. Os próprios Estados Unidos importam 
hoje mais 40% de petróleo do que em 1973, ano que marca o início da crise. 

Continuam a queimar carburante a um ritmo assustador, de tal forma que o consumo atingirá, este ano, mais 
de um bilião de toneladas. 

A falar com franqueza e contrariamente ao que se pensa, a crise da energia ainda não surgiu. Só a conhe- 
ceremos em 1985, e então com uma outra amplitude. Todos os especialistas mundiais concordam com esta previsão 
mais do que pessimista. 


1985 — O ano da supercrise 


É de resto fácil de compreender. O consumo mundial de petróleo é da ordem de um bilião e meio de tone- 
ladas. Ao ritmo actual, atingirá os dois biliões e meio de toneladas em 1985. E 1985 está à porta. 

Entretanto os países da OPEP (os exportadores do ouro negro) já fizeram saber que os seus Governos, nessa 
data, não tencionam ultrapassar uma produção total de 1800 milhões de toneladas. A crise é, portanto, inevitável. 

Em 1985, o consumo do mundo não comunista será asegurado 8% pela hidráulica, 18% pelo carvão, 20 % 
pelo gás natural e 8% somente pela energia nuclear. À fonte de energia essencial continuará, portanto, a ser petró- 
leo, com 46%, contra os 53 % actuais. 

Só daqui a alguns decénios é que o progresso técnico abrirá a porta de novas energias, termo-nucleares, 
mares, interior da terra e sol, 


Minicentrais solares móveis 


De todas as energias «paralelas», a fornecida pelo sol é a que tem mais possibilidades de se desenvolver 
rapidamente, o que, até agora, não era claro. É preciso reconhecer que as habitações solares (as casas aquecidas 
graças ao calor do astro) nada têm de estético. Quanto aos aquecimentos de água solares, a sua utilização é 
bastante limitada. 

Por outro lado, as centrais solares, previstas aqui e acolá, são verdadeiros monstros da técnica, com milha- 
res de hectares cobertos de painéis transmitindo o calor recebido para uma caldeira colocada no alto de uma torre 
de 150 metros. Que «beleza» na paisagem! Mas, eis que de repente, tudo vai mudar. Mini-centrais nucleares com 
potência de um ou dois kilowatts, apresentarão os seus espelhos curvos e reclináveis sobre uma circunferência de 
seis metros de diâmetro apenas. A caldeira, destinada a produzir o vapor, estará somente a dez metros de altura. 
Será o ideal para as aldeias do Terceiro Mundo ou para os bairros das nossas cidades. Fala-se mesmo já de mini- 
-centrais móveis colocadas sobre camiões. 

As descobertas vão em bom andamento. Na Austrália, na Unversidade de Sydney, para captar o calor solar 
com muito mais eficácia (a uma taxa de absorção de 809%) foi criada uma bateria de tubos de vidro com um reves- 
timento metálico, e colocado no vazio. Este processo será comercializado a partir de 1985. 

O sol não só dá calor, como também luz. É esta última que as pilhas solares transformam directamente 
em electricidade. Hoje reservam-se sobretudo para os satélites, porque custam extremamente caro (mais de 1500$00 
por watt). É que são fabricadas a partir de finíssimas lamelas de silício talhado num monocristal, cuja produção é 
muito onerosa. 

Dentro de 20 anos, porém, a coisa será diferente graças a um novo material, o silício amorfo que em cama- 
das muito finas, será colocado em suportes relativamente em conta. E o sol será cem vezes mais barato. 


(Informação do «Expresso», de 9-2-1980) 
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BORGES DA SILVA, J. F— O conceito de deslo- 
camento eléctrico na teoria do campo electroma- 
gnético de Maxwell, «Técnica», Lisboa, LIV (455) 
Dez. 1979, p. 133-140. 


Baseados no relato do próprio Maxwell no seu céle- 
bre Tratado da Electricidade e Magnetismo, faz-se uma 
tentativa de seguir o raciocínio físico que conduziu aos 
conceitos de Deslocamento Eléctrico e da Corrente de 
Deslocamento que lhe está associada. Estes conceitos 
desempenharam um papel decisivo na formulação das 
leis de Maxwell do campo electromagnético, que se 
mantêm ainda hoje, mais de um século depois da sua 
descoberta, como uma das grandes conquistas da fisica 
teórica. 

A circunstância de se verificar que algumas hipóte- 
ses chave, fundamentais para a génese da teoria, já 
não são sustentáveis à luz de conhecimentos mais re- 
centes, constitui um fenómeno interessante para os 
historiadores e filósofos da ciência ponderarem. 


ANDRADE E SILVA, João de-— 4 obra cientilica 
de Albert Einstein, «Técnica», Lisboa, LIV (455) 
Dez. 1979, p. 141-144, 


FIGANIER, João — Relato sucinto do Il Congresso 
Internacional! para o ensino da engenharia [Darms- 
tadt, Outubro de 1978), «Técnica», Lisboa, LIV 
(455) Dez. 1979, p. 145-150. 


FORTES, M. Amaral — Efeito de diversos parâme- 
tros nas propriedades medidas em ensaios de trac- 
ção, «Técnica», Lisboa, LIV (455) Dez. 1979, p. 
151-157. 


Com base numa teoria fenomenológica da deforma- 
ção plástica procura precisar-se o significado das pro- 
priedades mecânicas medidas em ensaios de tracção e 
estuda-se o efeito de diversos parâmetros — designada- 
mente a velocidade do ensaio, as dimensões do provete 
e a rigidez da máquina — nessas propriedades e, de um 
modo geral, nas curvas de tracção correspondentes à 
parte homogénea da deformação. 


E ge E A SE SS E 


SILVA JOÃO, A. F., MELO ARRUDA, A. M. e 
NIETO DE CASTRO, C. A.— Propriedades de 
transporte de gases diluídos, «Técnica», Lisboa, LIV 
(455) Dez. 1979, p. 159-167. 


Faz-se uma introdução à teoria de CHAPMAN-ENS- 
KOG para gases a baixa pressão, aplicando-a ao estudo 
da viscosidade, da condutibilidade térmica e do coefi- 
ciente de autodifusão de moléculas simples (Kr, Ns, 
CH) e ao coeficiente de difusão mútua duma mistura 
d Ar-Kr. Esta aplicação requereu o cálculo dos parâme- 
tros moleculares ; e q para os vários gases, inerentes 
a três modelos de potencial intermolecular, os modelos 
de Lennard-Jones 6-13, Dymond-Alder e de Kihara, o 
que foi efectuado pelo método das rectas isotérmicas 
e dos estdaos correspondentes, Faz-se uma discussão 
da validade dos métodos aplicados para a previsão das 
propriedades de transporte dos gases diluídos, compa- 
rando-se os parâmetros obtidos com os obtidos pelas 
técnicas de inversão e de ajuste a várias propriedades. 


MOURA, José Manuel Fonseca de — À teoria dos 
sistemas no contexto do controlo e da comunicação, 
«Técnica», Lisboa, LIV (455) Dez. 1979, p. 169- 
-182. 


Estuda-se neste relatório a disciplina Sistemas de 
Controlo e de Comunicação do grupo «Análise e Teoria 
dos Sistemas» do curso de Engenharia Electrotécnica. 
Começa-se por reflectir sobre o âmbito da Teoria dos 
Sistemas, especificando para o grupo de disciplinas as 
matérias afins e o espirito prevalecente, Percorre-se 
depois sumariamente aspectos do ensino da Teoria dos 
Sistemas no curso de Engenharia Electrotécnica no 
IST, nomeadamente no seu Departamento de Teoria 
dos Sistemas. Finalmente discute-se a natureza da dis- 
ciplina de Sistemas de Controlo e de Comunicação con- 
cretizando-se o respectivo programa, o conteúdo, mé- 
todos de ensino e aspectos bibliográficos. . 


CORTEZ, L.. MUGE, F. e PEREIRA, H. — Planea- 
mento mineiro de um depósito de fosfatos, «Téc- 
nica», Lisboa, LIV (455) Dez. 1979, p. 183-187. 


Nos estágios iniciais do planeamento de uma nova 
exploração mineira, o erro cometido na estimação das 
reservas é um factor importante que não deve ser des- 
prezado. O cálculo deste erro por métodos Geoestatis- 
ticos permite a execução de uma análise de risco 
fixando certas condições econômicas, 

É descrito no artigo um exemplo concreto de apli- 
cação destas técnicas. Trata-se de um depósito de fos- 
fato explorado a céu aberto para o qual foram obtidas 
curvas tonelagem/teor com os respectivos intervalos de 
confiança. Estabeleceu-se uma relação entre estas cur- 
vas e uma função económica para três malhas de reco- 
nhecimento, tendo as subsequentes decisões sido toma- 
das com base num balanço entre o custo e o valor da 
informação suplementar. 
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SILVA JOÃO, A. F., MELO ARRUDA, A. M. 6 
NIETO DE CASTRO, C. A. — Transport properties 
of dilute gases, «Técnica», Lisboa, LIV (455) Dez. 
1979, p. 159-167, 


This paper presents an introduction to the CHAP- 
MAN-ENSKOG theory of dilute gases, and an applica- 
tion of it to the study of viscosity, thermal conducti- 
vity and self diffusion of simple molecules (Kr, Na, 
CH and to the mutual diffusion coefficient of the 
Ar-Kr mixture. The estimation of the intermolecular 
potential constants & and has been made for the 
several gases, using the Lennard-Jones, Dymond-Alder 
and Kihara models, and two methods, the isotherms 
method and the corresponding states principle one. 
A discussion of the application of this theory to the 
prevision of transport properties of dilute gases is also 
made, followed by a comparison of the molecular para- 
methers obtained with those obtained by inversion 
techniques of fitting procedures. 


MOURA, José Manuel Fonseca de — Systems 
theory in the context of control and communication, 
«Técnica», Lisboa, LIV (455) Dez. 1979, p. 165- 
-182. 


The paper presents Control and Communication 
Systems, a subject of the group «Analvsis and Theory 
of Systems» of the Electrical Engineering program. If 
starts out by reflecting on the nature of Systems 
Theory itself, specifying for the above group of sub- 
jects its contents and prevalling attitude, The report 
then summarizes aspects of the teaching of Systems 
Theory in the Elecirical Engineering Course at IST, 
namely at its Systems Theory Department. Finally, 
the paper discusses Control and Communication Sys- 
tems presenting its programm, its contents, the peda- 
gegical methods used and the required bibliography. 


CORTEZ, L. MUGE, F. e PEREIRA, H. — Mine 
planning in a phosphate deposit, «Técnica», Lisboa, 
LIV (455) Dez. 1979, p. 183-187. 


In the early stages of the planning of a new 
exploitation, the error in ore reserve estimation is an 
important factor that must not be disregarded. As this 
error can be computed by geostatistical methods, the 
risk of the venture can be calculated for different 
recognition patterns in a certain asumption of econo- 
mical conditions. 

A case story is described for an open pit phosphate 
deposit in which grade/tonnage curves and their confi- 
dence intervals are related to an economic function for 
three drilling patterns and decisions are taken on the 
basis of cost and value of supplementary information, 
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BORGES DA SILVA, J. F.——- The concept of elec- 
tric displacement in Maxwell's electromagnetic field 
theory, «Técnica», Lisboa, LIV (455) Dez. 1979, 
p. 133-140, 


Based on Maxwell's own account in his celebrated 
Treatize on Electricity and Magnetism, an attempt is 
made to follow the physical reasoning leading to the 
concepts of Electric Displacement and its associated 
Displacement Current, These concepts played a decisive 
role in the formulation of Maxwell electromagnetic field 
laws, which still stand today, more than a century 
after their discovery, as one of the greatest achieve- 
ments of theoretical physics, 

The circunstance of it being found that some key 
assumptions, fundamental to the genesis of the theory, 
are no longer tenable in the light of more recent 
knowledge, is an interesting phenomenon for historians 
and philosophers of science to ponder, 


ANDRADE E SILVA, João de — Cientifico work of 
Albert Einstein, «Técnica», Lisboa, LIV (455) Dez. 
1979, p. 141-144. 


FIGANIER, João — A brief relate of the Il Interna- 
tional Congress on Engineering Teaching (Darms- 
tadt, October 1978), «Técnica», Lisboa, LIV (455) 
Dez. 1979, p. 145-150., 


FORTES, M. Amaral — Effect of various parameters 
on the properties measured in tensile tests, «Téc- 
nica», Lisboa, LIV (455) Dez. 1979, p. 151-157. 


Using a phenomenological theory of plastic defor- 
mation an analysis is made of the significance of the 
various mechanical properties obtained from tensile 
tests and of the effect of parameters, such as the strain 
rate, the dimensions of the specimens and the tensile 
machine stifness, on these properties and, in general, 
em pç tensile curves for the homogeneous deformation 
n tension. 
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